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ABSTRACT 

In this paper we estimate binary compact object merger detection rates for LIGO, including the 
potentially significant contribution from binaries that are produced in elliptical galaxies near the epoch 
of peak star formation. Specifically, we convolve hundreds of model realizations of elliptical- and spiral- 
galaxy population syntheses with a model for elliptical- and spiral-galaxy star formation history as a 
function of redshift. Our results favor local merger rate densities of 4 x 10~^ Mpc~^ Myr~^ for binary 
black holes (BH), 3 x 10"^ Mpc"^ Myr~^ for binary neutron stars (NS), and 10"^ Mpc"^ Myr"^ for 
BH-NS binaries. We find that mergers in elliptical galaxies are a significant fraction of our total 
estimate for BH-BH and BH-NS detection rates; NS-NS detection rates are likely dominated by 
the contribution from spiral galaxies. Limiting attention to elliptical- plus only those spiral-galaxy 
models that reproduce current observations of Galactic NS-NS, we find slightly higher rates for NS-NS 
and largely similar ranges for BH-NS and BH-BH binaries. Assuming a detection signal-to-noise ratio 
threshold of 8 for a single detector (in practice as part of a network, to reduce its noise) , corresponding 
to radii Dbns of the effective volume inside of which a single LIGO detector could observe the inspiral 
of two IAMq neutron stars of 14Mpc and 197Mpc, for initial and advanced LIGO, we find event rates 
of any merger type of 2.9 x 10~^ - 0.46 and 25 — 400 per year (at 90% confidence level), respectively. 
We also find that the probability Pdctcct of detecting one or more mergers with this single detector can 
be approximated by (i) Pdotoct — 0.4 -I- 0.5 log(r/0.01 yr), assuming Dbns = 197 Mpc and it operates 
for T years, for T between 2 days and 0.1 yr; or by (u) Pdctcct — 0.5 + 1.5 log(Dbns/32 Mpc), for one 
year of operation and for -Dbns between 20 and 70 Mpc. 

Subject headings: Stars: Binaries: Close; Stars: Pulsars: General; Gravitational waves 



1. INTRODUCTION 

Over the last decade the question of the Galactic 
inspiral rate of binaries with two compact objects (neu- 
tron stars NS or black holes BH) has attracted attention 
primarily because of the development and planning of 
gravitational-wave interferometric detectors both on the 
Eart h and in space (e.g., LIGO and GEO600, described 
in lAbbott et al. (The LIGO Scientific CoUaborationil 
(j2003| ): Virgo, described at the Virgo project website 
www.virgo.infn.it; and LISA, at lisa.nasa.gov). 
These rate estimates have been widely used in the 
assessment of gravitational inspiral detectability, 
given assumed instrument sensitivities. A number 
of different groups have calculated inspiral rates 
using population synthesis calculat ions, most com- 
monly with Monte Carlo m et hods (iFrver et al.l 119981: 
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iVoss fc TaurisI 120031 ). Such studies consider the 
complete formation history of double compact objects 
through long sequences of binary evolution phases, 
terminated by a gravitational-wave driven inspiral 
towards merger. Though less critical for studies of the 
Milky Way, contributions to the present-day merger 
rate per unit volume from early star formation has most 
often been ignored in merger rate calculations. Since 
long inspiral delays before merger are not sufficiently 
uncommon, even for double neutron stars, merger de- 
tection rate calculations should account for the full time 
history of star formation, inclu ding star formation in the 
early universe; see for example Ide Freitas Pacheco et al.l 
()2006D . Additionally, our current understanding of 
single and binary star evolution is incomplete. The 
many uncertaintites have been parameterized and the 
resulting parameter space explored in some studies 
to determine the range of plausible results (see, e.g., 
IBelczvnski et al.l l2008b[ ). both by focusing narrowly 
on physically- motiv a ted reg ions of parameter space as 
in IBelczvnski et al.l ()2008b[) and by broadly exploring 
all plausible population synthesis simulations (see, e.g., 
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lO'Shaughnessv et alj l2008cJ . and references therein; 
henceforth denoted PS-GRB). Finally, by assigning 
equal prior likelihood to any population synthesis 
model, PS-GRB estimated the relative likelihood of any 
BH-NS or NS-NS merger ratcQ 

The population of binary black holes, however, behaves 
qualitatively differently from BH-NS or NS-NS binaries: 
(i) most of those that merge do so only after a few to sev- 
eral Gyr delay after their birth as binary stars; and (ii) 
those BH-BH binaries with merger times from 1 — 10 Gyr 
have masses quite different from the typical assumption 
of a IOMq BH; see Appendix. The rarity, long delay 
times, high masses and thus cosmologically significant 
detection horizons of binary black hol es have already 
been discussed in the literature (see.e.g-. lKulczvcki et alJ 
llOOG.) and references therein). However, population syn- 
thesis simulations extensive enough to contain a statisti- 
cally meaningful binary black hole population come at a 
significant computational cost. Even with adequate pop- 
ulation synthesis simulations, given the long-lived nature 
of the progenitors, BH-BH merger and detection rates 
depend critically on the earliest and least certain star 
formation rates. Certain exotic but not exceptionally un- 
common binary evolution scenarios could lead to a "high 
rate tail" ~ a scenario with rare massive mergers, for ex- 
ample, that could plausibly produce a merger detection 
in the near future. In this paper we strive to understand 
what realistic scenarios could lead to a "high rate tail" 
that LIGO could detect or, conversely, constrain. 

Given the technical challenges and significant uncer- 
taint ies involved, only one paper (jO'Shaughnessv et al.l 
|2005[ ) has previously estimated the relative likelihood of 
different BH-BH detection rate estimates expected from 
isolated binary star evolution, only for the Milky Way. 
In the present study we calculate the expected detec- 
tion rate for LIGO while for the first time simultane- 
ously including (i) all past star formation, particularly 
the overwhelming importance of elliptical galaxies; (ii) a 
large database of results that accounts for the dominant 
model uncertainties; and (iii) a careful treatment of the 
mass distribution of BH-BH binaries. Additionally, un- 
like previous analyses presenting distributions of merger 
rates, we have post-facto varied the birth binary fraction 
of stars from 100% to 15%; this addition allows us to 
better compare our results with the fiducial population 
synthesis model and BH -BH merger rates presented in 
iBelczvns ki et al.l ()2007al ). Finally, whereas previous com- 
parisons relie d essentially on confidence in tervals (e.g., 
PS-GRB and lO'Shaughnessv et al.l (|2008bD . henceforth 
denoted PSC2), in this paper we generate true posterior 
distributions: each model is weighted by its relative con- 
ditional likelihood given observations of merging NS-NS 
binaries in the Milky Way. 

To fully assess the total probability for a LIGO detec- 
tion, rather than limit attention to BH-BH mergers we 
explore the rate at which all LIGO-detectable binaries 
(BH-BH, BH-NS, and NS-NS; henceforth collectively de- 

1 Previously, lO'Shaughnessv et al.l 120051. I2008bl ') estimated 
LIGO binary merger detection rates by extrapolating from the local 
Milky Way based on the blue light density in the nearby Universe. 
As we discuss here, their estimate is accurate for short distances 
and for binary merger progenitors that are not preferentially long- 
lived (i.e.. the median delay betwreen birth and merger is less than 
IGyr). 



noted as compact binary coalescences or CBCs) merge 
through gravitational wave emission. As the methods 
used to determine LIGO detection rates from population 
synthesis simulations a re already extensiv ely discussed in 
the literature fsee.e.g. iBulik et al .112001 and PS-GRB), 
our presentation only reviews those tools, emphasizing 
unique features of our present analysis (i.e., our methods 
to estimate BH-BH masses). Specifically, in §[2] we briefly 
review the ingredients that enter into an estimate of the 
gravitational wave detection rate for a short-ranged net- 
work {zmax ^ !)• While the local universe is emphasized 
in the text, networks of advanced ground-based interfer- 
ometers can detect optimally oriented BH-BH binaries 
with component masses M ~ 10 — ISA/q at cosmologi- 
cally significant distances z > 0.1. Merger rates on the 
past light cone of a detector become inhomogcneous (ver- 
sus redshift) at this scale, simply because the star forma- 
tion rate increases dramatically near z ~ 1 — 2 during the 
epoch of galaxy assembly. At these distances, the detec- 
tion rate must be integrated over the mass distribution, 
the full networked orientation-dependent sensitivity, and 
redshift-dependent merger rate. This integral becomes 
increasingly sensitive to rare, high-mass events that can 
be seen in the ancient universe. Therefore, we general- 
ize the short-range method described in § [5] to include 
all sources on the past light cone, to arbitrary redshift 
and with accurate orientation-dependent sensitivity. In 
§ [3] we review the model for star formation in the uni- 
verse adopted here. This experience is applied in §|4](for 
BH-BH binaries) and § [5] (including NS-NS and BH-NS 
binaries too), where we review our population synthesis 
calculations; and explain what features of those calcula- 
tions influence our predictions for the relative likelihood 
of a LIGO detection. These sections also describe Fig- 
ure [71 which contains our predictions regarding initial 
and advanced LIGO detection rates. Finally, in § 15.21 
we show how rate constraints from the observed sample 
of Milky Way NS-NS binaries affect the predictions for 
LIGO detection rates. 

To summarize, by fully simulating the past history 
of the local universe, this paper develops models for 
the present-day detection rate of short-range (z ^ 1) 
and long-range {z not much less than 1) gravitational- 
wave detectors. Our results are rate distributions, where 
each distribution includes some normalization uncertain- 
ties (star formation rate and fraction of stars born in 
binaries); certain population synthesis model parame- 
ters; and our simulation Monte Carlo uncertainty. Fur- 
ther, our Bayesian approach to model constraints surveys 
many key uncertainties that should be included, what- 
ever the model family involved, when attempting to in- 
terpret upper limits or detections from pulsar popula- 
tions and gravitational wave observatories. We believe 
our distributions represent the best predictions of a con- 
crete, conservative model family that includes both ellip- 
tical and spiral star formation yet is also consistent with 
initial LIGO upper limits on CBC merger rates. 

2. GRAVITATIONAL WAVE DETECTION 

This section primarily reviews the final stage in 
any calculation of a LIGO detection rate: the con- 
nection between, on the one hand, the event rate per 
comoving volume TZ and mass distribution p(mi,m2) 
of merging binaries and, on the other hand, the de- 
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tection rate Rjj [Eq. ©]. This relation depends 
critically on the range to which LIGO could observe 
each merger. A thorough discussion of the LIGO 
range requires careful review of data analysis strategies 
and interferometer network geometry, and remains 
subst antially beyond the scope of this paper. For details 



see Abbott et al. (The LIGO Scientific CollaborationH 



200£ ): IAbbott et al fThe LIGO Scientific Collaboration)! 



20091) and references therein. For example, with mod- 



ern multi-mass-region searches using inspiral templates, 
one essential effect is a merger detection threshold 
in SNR that changes with total binary mass, as the 
well-understood portions of detection templates grow 
shorter. In this paper, we adopt the customary set of 
approximations (described below) and ignore the errors 
they introduce. 

The distance to which LIGO is sensitive can be calcu- 
lated using the sensitivity of the detector (strictly, the 
strain noise spectral distribution Sh), the masses of the 
inspiralling bodies, and the emitted waveform. [For sim- 
plicity, we will assume the black holes initially nonspin- 
ning.] Specifically, for a source with masses nii and 7TI2, 
merging at a luminosity distance D (or equivalently at 
redshift z, or comoving distance r) with a binary incli- 
nation i and a relative orientation to the detector given 
by 9, (j) (the orientation on the sky) and ip (a phase angle 
in the plane of the detector), the signal-to- noise seen by 
a matched-filter search in the LIGO data stream can be 
expressed relative to the strain noise spectrum of a sin- 
gle interferometer Sh and the Fourier transform h{f) of 
gravitational wave strain h{t) received at the detector by 
an arbitrarily oriented and located source by 



|rfe(mi(l + z), 7712(1 + z),/)p 
Sh{f) 

(1) 

using units with G = c = 1, the known dependence 
of planar-symmetry gravitational waves on angle, and 
a standard one-sided Fourier convention for h and Su 
versus frequency /. The function w{6, (pjip, l) takes val- 
ues between and 1 and completely encompasses the 
detector- and source-orientation dependent sensitivity 
for the most common sources, those that are dominated 
hy I = m = 2 quadrupole emissior0 

wi9, (/), V-, l)^^FI{1 + cos2 /,)2/4 + F2 cos2 t (2) 



= -(1 + cos 9) cos 2(/) cos 27/- 
+ cos 9 sin 20 sin 2ij} 
Fx =-i(l + cos^ 9) sin 2(/) cos 2^ 



(3) 



-cos^?cos2(/)sin2'0 



(4) 

To a rough approximation that depends on the data 
analysis strategy used and the amount of nongaus- 
sian noise present in the detector, a single LIGO 
interferometer can detect the gravitational wave sig- 
nat u re of a merging binary if p > Pc (see, 
e.g., lAbbott et al. (The LIGO Scientific Collaboration)! 

^ Peak sensitivity is attained for a source directly overhead with 
orbital plane normal to the line of sight. 



()2006l[2008l ). and references therein; henceforth we adopt 
Pc = 8). For this reason, the previous expression for the 
signal to noise [Eq. ([ij] is often re-expressed as 



p = PcW{9,(l),'lp,L 



£'H(mi(l + z), 7772(1 -fz)) 
D 



(5) 



which, by comparison with Eq. ([T]), implicitly defines the 
horizon distance Dn - the maximum luminosity distance 
to which the detectors are sensitive - as a function of the 
redshifted mas ses 77-71 (1 + z) and 7712 ( 1 -I- z) of the bi- 
nary ( see, e.g.. [Anderson et all (12001 ). iFinn fc Chernofi! 
(fT99l and [Flanagan fc Hughed fTool for a brief review 
of the theory underlying the direction-dependent LIGO 
sensitivity as well as for expressions that approximate 
Dui'mi,'m2) for low-mass CBC binarief0). 

More commonly used is the volume-averaged distance 
Dy , chosen so the volume of a sphere of radius agrees 
with the average volume enclosed in a (Euclidean) detec- 
tion surface: 



1 

An 

:(i?H/2.26)3 



,^ d'ib d cos L o 
dQ — w'^ 



(6) 



The volume averaged distance is a meaningful measure of 
sensitivity only when the horizon range is much smaller 
than the Hubble scale. 

A network of detectors can coherently add signals from 
each interferometer to increase the signal-to-noise p asso- 
ciated with each signal and by implication its reach. To a 
first approximation, ignoring small differences in sensitiv- 
ity due to both their intrinsic differences and their orien- 
tations|3 the combined SNR all interferometers is higher 
by roughly \/J2k ^k/^i where Lk is the length of the 
fcth interferometer. For example, this factor is approxi- 
mately gN ^1 + 1 + 1/4 ~ V2.25 for the initial LIGO 
network , consisting of two 4 km ari d one 2 km interferom- 
eter; see lCutler fc Flanaganl (|1994( ) for details on realistic 
multidetector beampatterns.. However, in part because 
a network performs more trials of the same data streams 
(e.g., for each sky position) and is sensitive to both po- 
larizations simultaneously, for the same false alarm rate 
the network detection threshold Pc.net is generally greater 
than the threshold Pc for the indivi dual detectors; see for 
examp le ICutler fc Flanagajil (|1994[ ) and IBuonanno et al.! 
(|2003f ) for a discussion in the case of gaussian noise. 
One can therefore speak of a single-interferometer (Z3h) 

and network (£'h''"' — 9 nDyi[Pc/ Pc.net)) horizon dis- 
tance. Unfortunately, only expert analysis of real data 
can determine the sensitivity of real gravitational-wave 
networks, not the least because noise statistics (and 

^ A more thorough review of the orbital geometry and waveforms 
underlying the reach of a single-interferomet er s earch for circularly 
inspiralling point particles can be found in Brown ( 2 0041) and (for 



the more c omplex case of spinn ing binaries) in |A£ostolatos_e^^ 
(jT994); Bu onanno et al.l lj2003 ). Additionallv. IKopparamT'et'ld 
112008, ) discusses how this orientation-dependent sensitivity influ- 
ences the present-day LIGO search. 

* The small bias (~ 10%) introduced by this ap- 
proximation is much smaller than the characteristic un- 
certainties discussed later in this paper; see, for exam- 
ple, LIGO site locat ion inf ormation in LIGO-T980044- 
10 (available at |http: / / admdbsrv.ligo.caltcch.cdu/dcc/1 
and the LAL software do cumentation ( avai lable tr om 
|http: / /www.lsc-group.phys.uwm.edu/daswg/projccts/lal.html[ |. 
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therefore Pc,net) are highly detector- and even search- 
dependent. To avoid ambiguity and misleading approxi- 
mations, we provide results for single, idealized intcrfcr- 
omctcrs assuming a detection threshold of = 8 (e.g., 
appropriate to gaussian noise). This reference sensitivity 
agrees with the customary sensitivity measure used in 
the LIGO coalescence search, a threshold of pc = 8 for a 
single detector as part of a network For sources in the 
local universe, without loss of generality the reader can 
scale our results to any realistic network and search. 

For sources of sufficiently low mass {Mtot ^5 20Mq), 
LIGO is primarily sensitive to less-relativistic and there- 
fore better-understood phases of the spiral-in, allow- 
ing u s to approxim ate the waveform by its Newtonian 
limit (|Peterslll96^. After considerable alg ebra (see, e.g. 
iFlanagan fc HughedlTOOl ICutler fc Flamg an 199|), the 
signal-to-noise due to an optimally-oriented inspiral at 
comoving distance r (corresponding to a redshift z{r)) 
can be related to the strain sensitivy Sh by 



P 



df 



Shif) 



(7) 



where A4c = (mi 7712)'^/''' /(mi + 1122)^^^ is the "chirp 
mass" of the binary. Based on this expression, individ- 
ual interferometers in the initial , and advancecQ LIGO 
networks can see low-mass binaries (toi +m2 ^ 4OAfQ)0 
out to a horizon (volume-averaged) distance 



7^H(v) =^Ch(v)(A^c/1.2A/0)^/^ 



(8) 



where Ch^v) = 31(14) Mpc and Ch(v) = 445(197) Mpc 
respectively. For clarity, the distance Cy will henceforth 
be denoted Dbns, the volume-averaged distance to which 
a binary neutron star inspiral can be detected by a single 
interferometer as part of a network. 

In terms of this horizon distance Dn, the intrinsic 
present-day rate of events per unit comoving volume 
7^.(0), and the chirp mass distribution of merging binaries 
p(A^c), the average rate of events with p > pc occurring 

For the real LIGO interferometers, no single detector search at 
this threshold is possible; multi-detector coincidence is required to 
reduce the nongaussian background to roughly the level indicated 
(hence "as part of a network" ) . 

® We use a published LIGO S5 sensitivity (for the 
LHO interferometer, on March 13th, 2006), available at 
|http: / / www . ligo . caltech . edu/ ~ j z weizig / distr ibution/LSC_Data[ 

' Wc adopt an advanced LIGO noise curve Irom 
LIGO T0900288, available as ZERO_DET_high_P.txt at 
[dccdigo.org/DocDB/0002/T0900288/002/ This noise curve 
is taken from the GWilNCJ program with Mode lb parameters, as 
described in LIGO document LIGO-T070247-01-I: 125 W input 
laser power, 20% signal recycling mirror (SRM) transmissivity, 
and no detuning of the signal recycling cavity. 

* For comparison, full strong-field numerical relativity waveforms 
suggest a single-interferometer initial LIGO horizon distance Dh 
to the most massive equal-mass binaries M ~ 40 that differs from 
the low-mass limi t by 2%, multiplied by a factor of order unity; see 
lAiith et al.l II2009I ). On the contrary, our limitation to equal-mass, 
nonspinning binaries introduces greater error: for binaries M < 
I6M0, mass ratio corrections are roughly 4%, while black hole 
spin-dependent corrections are O.lJ/m^^, multiplied by a factor 
of order unity (if aligned). At the extreme, BH-NS binaries with 
strong spin-orbit misalignment and thus precession could be visible 
only out to O.SDh^ej, multiplied by a factor of order unity, for 
Duref the horizon distance to a comparable aligned binary. Such 
extreme uncertainties in range are less plausible, as isolated binary 
evolution models favor spin-orbit alignment. 



in the nearby universe can be expressed as 

-p{Mc)dMc (9) 



RD=n{0) f 



9 , ,r^dib d cos L 
r^drdfl—- 



7^(0) 



J dMcYDv{Mcfp{Mc) 



(10) 



where the integral is taken over all detectable combina- 
tions of masses, orientation, and location. For low-mass 
binary mergers {M < 40A'/o), for which Eq. ^ applies, 
the chirp mass-weighted average simplifies to 

i?z5^7^(0)^^^ J p{Mc)iMc/l.2MQy'^'dMHl) 

an expression that has been applied extensively in 
the literature to translate event r ates per unit volume 
into L IGO detection rates (see, e.g. lO'Shaughnessv et al.l 
l2008bl and references therein). 

2.1. Rates for advanced detectors 

As the product of the present-day merger rate and an 
effective detection volume, the detection rate estimate 
Eq. ^ assumes a locally homogeneous and isotropic uni- 
verse. Advanced detectors can see back to sky-position- 
dependent epochs with noticeably different star forma- 
tion history {z > 0.1). For advanced detectors, a more 
generic expression must be used, which integrates over 
the mass distribution, the redshift-dependent merger 
rate, and the orientation-dependent reach w. 

For the purposes of this paper, we continue to assume 
the maximum luminosity distance to which a single inter- 
ferometer can identify a merger in noise is well-described 
by the inspiral phase [Eq. [7] . To normalize our estimate, 
we adopt a singlc-intcrferomcter advanced LIGO range 
of 445 Mpc. Based on the optimal range for a given chirp 
mass, each simulation's the chirp mass distributions, and 
the single- interferometer beampattern w [Eq. [2], we de- 
termine the fraction of all mergers at redshift z that could 
be detected, Pok{z): 



Pok{z)-^ 



dfl dtp d cos t 

D{z)<wDn TT 2 



dMcPiMc) (12) 



The overall single-interferometer detection rate is there- 
fore the sum over the past light cone of the (redshifted) 
rate of mergers on it, times the fraction Pok of mergers 
that could be detected: 



Rd= / dz 



dv nit) 

rfz 1 + z 



Pokiz) 



(13) 



For a maximum reach z <C 1, these two expressions are 
equivalent to Eq. [T] 

3. MULTICOMPONENT STAR FORMATION 
HISTORY 

The LIGO detection rate depends on the mass distri- 
bution of merging binaries p{Aic) and on the average 
rate of mergers per unit comoving volume, TZ{t). This 
rate, in turn, generally depends on the net contribution 
from multiple star-forming populations, a contribution 
that convolves the star formation rate in each component 
(the subject of the present section) with the rate at which 
each components' star-forming gas yields CBC mergers 
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Fig. 1. — Star formation history of the universe used in this 
paper versus time, relative to the pres ent day. Solid line: Net star 
formation history implied by Eq. I I16II . Dashed, dotted line: The 
star formation history due to elliptical and spiral components. 

(the subject of §[3]). More specifically, the merger rate 
density TZ{t), a quantity required to calculate the detec- 
tion rate, is obtained from (i) the star formation rate in 
each component dpc/dt where C indexes the different 
star-forming types, C — e for ellipticals and s for spi- 
rals, defined as the mass per unit comoving volume and 
time that forms as stars; (ii) the mass efficiency Ac at 
which each type of CBC binaries form, defined as the to- 
tal number of binaries that survive isolated evolution and 
form CBC binaries per unit star forming mass; and (iii) 
the probability distribution dPmc/dt for merger events 
to occur after a delay time t after star formation, where 
dPmC is the fraction of mergers occuring between t and 
t + dt after the progenitor binary forms at i = 0: 

7^(^) = ^7^c(^) (m) 

c 

licit) = j'^dhXc^{t^t,)^{t,) (15) 

where the integration variable tt is time at which Xdp 
bi naries per unit volume are born . As first recognized 
bv lde Freitas Pacheco et ahl (|2006[ ) and as demonstrated 
systematically in PS-GRB, a multicomponent star for- 
mation history that includes both elliptical and spiral 
galaxies must be applied when calculating present-day 
merger rates of double compact objects. Simply put, 
even though elliptical galaxies formed stars long ago, 
dP/dt decays slowly enough (roughly as 1/t; see, e.g., 
PS-GRB) that the high rate of star formation in the 
early universe can lead to a significant clliptical-galaxy- 
hosted merger rate density in the present-day universe. 
Because elliptical galaxies have distinctly different star 
forming conditions (e.g., metallicities) than present-day 
star-forming galaxies, models for stellar evolution in the 
universe should not require identical behavior in present- 
day solar-metallicity galaxies and in young ellipticals. A 
full multicomponent treatment is necessary. The first 
multicomponent star formation history was applied to 
predict BH-NS and NS-NS merger rates in PS-GRB; this 
paper employs the same framework and preferred star 
formation history, discussed below. 
We adopt the two-component star formation history 



model presented bv lNagamine et ahl (|2006[) . This model 
consists of an early "elliptical" component and a fairly 
steady "spiral" component, with star formation rates 
given by 

p = Pe + Ps (16) 

pc = Acit/Tc)e-'/^^ (17) 

where cosmological time t is measured starting from 
the beginning of the universe and where the two 
components decay timescales are re,s = 1.5 and 
4.5 Gyr, respectively (s ee Section 2 and Table 2 
of iNagamine et al.l |2006[) . These normalization con- 
stants A^.,, = . 22,0.0 6M(T^yr~^Mpc~^ were chosen by 
INagamine et al.l ()2006[ ) so the integrated amount of el- 
liptical and spiral star formation reproduce the present- 
day census of baryonic matter in ellipticals and spirals, 
allowing for a certain fraction of gas recycling between 
different generations of starsQ Figure [T] illustrates the 
star formation rates assumed by the model adopted here. 

Each component forms stars in its own distinctive con- 
ditions, set by comparison with observations of the Milky 
Way and elliptical galaxies. We assume mass converted 
into stars in the fairly steady "spiral" component with 
solar mctallicity and with a fixed high-mass initial mass 
function (IMF) powe r law \p = —2.7 in the brok en-power- 
law Kroupa IMF; see lKroupa fc Weidnei] (|2003[ )] . On the 
other hand, we assume stars born in the "elliptical" com- 
ponent are drawn from a broken power law IMF with 
high-mass index within p E [—2.27, —2.06] and mctal- 
hcity Z within 0.56 < Z/Zq < 1.5. These elliptical 
birth conditions agree with obse rvations of bot h old el- 
lipticals in the local universe (see iLi et al.ll2006l and ref- 
cren ces therein) as well as of young sta rburst clusters 
fsee lFaU et al|[2005t [Zhang fc Fall|[l99l and references 
therein). 

Since two independent stellar populations give birth to 
CBCs, the net detection rate Rd is the sum of the de- 
tection rate due to ellipticals and spirals. In particular, 
the probability distribution function for Rjj is necessarily 
the convolution of the distribution functions for detection 
rates due to elliptical and spiral galaxies, respectively. 
Because of the many orders of magnitude uncertainty in 
merger rates and the need to smooth across relative er- 
rors, however, all our PDFs are stored and smoothed in 
a logarithmic representation, as p{\ogR]j); see for exam- 
ple PS-GRB. For this reason, the convolution takes on an 
unusual form, involving an integral over the fraction / of 
the detection rate due to ellipticals (i.e., / = RD,e/ Rd)- 
Using this variable to form the logarithmic convolution, 
the distribution for the net rate logi?D can be calculated 
from the distributions for ellipticals and spirals {jpe,s) as 

^^'-^''-^=1 /(l-.mn(10) 
xpe(log Rd +log/) 
xp,(logi?z5+log(l-/)) (18) 

= Ps®Pe (19) 

Limitations of two- component model: Our two compo- 
nent model implicitly adopts two extremely strong as- 

^ Our previous paper PS-GRB incorrectly listed a larger value 
As = O.15M0yr-iMpc-3. 
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sumptions: (i) that metallicity in each component does 
not significantly evolve with redshift and (ii) that all 
star-forming gas in each component is identical. While 
sufficient for most problems involving binary evolution, 
these approximations may be particularly ill-suited to 
determining the progenitors of gravitational wave detec- 
tions. Very low metallicity environments could produce 
BHs with exceptionally high mass; also, the BH binary 
merger rate increases strongly with metallicity (see, e.g. 
iBelczvnski et al.ll2009| ). Because of the long typical de- 
lay between progenitor birth and BH-BH merger, very 
low-metallicity environments in the early universe could 
contribute significan tly to the present - day gra vitational- 
wave detection rate IBelczvnski et al.l (|2008aD . Even at 
present, low-metallicity environments exist often enough 
to potentially dominate the local merger rate (see, e.g., 
lO'Shaughnessv et al.l 120083 ). For example, the nearby 
BH-BH progenitor binary IC 10 X-1 both lies in a low 
metallicity environment and suggests a high BH-BH de- 
tection rate for initial LIGO (0.5 yr, multiplied by a fac- 
tor of orde r unity, strongly de pendent on survey selection 
effects; see lBulik et all pOOl ). 

Despite strong theoretical and observational evidence 
for a significant low-metallicity contribution, we have 
neither the confidence in the necessary astrophysical in- 
puts (i.e., a redshift-dependent metallicity distribution) 
nor the required databse of simulations to include these 
metallicity effects at this time. Further, unless initial 
LIGO detects BH-BH mergers, the peak detection rates 
cannot be significantly larger than what our model fam- 
ily predicts; see Figure [T] For this reason we continue to 
adopt the simple and conservative two-component mech- 
anism outlined above. 

Relation to blue light density: Previously, merger rate 
estimates for the Milky Way have been scaled up 
to the local universe, assuming me rgers track blue 
light (see, e.g. iKopparapu et al.l 120081 and references 
therein). As rou ghly speaking blue light traces star for- 
mation Csee.e.g. iGa llagher et al.|[T984l : iKennicuttI 119981 : 
lO'Shaughnessv et a l. 2008c), this extrapolation roughly 
agrees with our estimate for binary merger rates in 
spiral-like galaxies, with a suitable choice of normal- 
ization (i.e., star formation rate per unit blue light); 
see t he conclusions for a detailed numerical compari- 
son. lO'Shaughnessv et al.l (|2008c|) explains why blue 
light normalization was not adopted in this paper. 

4. GBC POPULATION SYNTHESIS MODELS 

As outlined in §[2]and §[2][see particularly Eqs. (|9ll5p ]. 
in order to estimate how often CBC binary mergers 
could be seen with gravitational- wave detectors we must 
know several features of the CBC population: (i) how 
many CBC binaries form from any given star forming 
mass (characterized by the mass efficiency A); (ii) how 
long these CBC binaries last between birth as stars and 
merger through gravitational wave emission (character- 
ized by dP/dt); and (iii) what masses these merging CBC 
binaries have (characterized by a probability distribution 
in chirp mass p{M.c))- In the absence of observations, 
these properties must be obtained from theoretical mod- 
els of stellar and binary evolution. We study the for- 
mation of compact objects with the S tarTrack popula- 
tion synthesis code, first developed by IBelczvnski et al.l 
((2002i) and recently significantly extended as described 



in detail in IBelczvnski et al.l ()2008bD . Since our under- 
standing of the evolution of single and binary stars is 
incomplete, this code parameterizes several critical phys- 
ical processes with a great many parameters (~ 30), 
many of which influence compact-object formation dra- 
matically. In this specific study, in addition to the IMF 
and metallicity (which vary depending on whether a bi- 
nary is born in an elliptical or spiral galaxy), seven pa- 
rameters strongly influence compact object merger rates: 
the supernova kick distribution (modeled as the super- 
position of two independent Maxwellians, using three pa- 
rameters: one parameter for the probability of drawing 
from each Maxwellian, and one to characterize the dis- 
persion of each Maxwellian), the stellar wind strength, 
the common-envelope energy transfer efSciency, the frac- 
tion of angular momentum lost to infinity in phases of 
non-conservative mass transfer, and the relative distri- 
bution of masses in the binary. For this reason, fol- 
lowing methods largely outlined in our previous work 
in PS-GRB and PSC2, we perform a broad parameter 
study to ensure all plausible population synthesis mod- 
els have been explored. We include simulations of two 
different classes of star-forming conditions: "spiral" con- 
ditions, with Z = Zq and a high-mass IMF slope of 
p = —2.7, and "elliptical" conditions, with a much flat- 
ter IMF slope p ~ —2.3 and a range of allowed metal- 
licities 0.56 < Z/Zq < 1.5. Though we have performed 
thousands of simulations of these two types of conditions 
in the past (e.g., PS-GRB), to obtain better statistics on 
rare BH-BH events, we have constructed a new set of 
archives of Ns = 282 and Ne = 206 simulations of 
these two type of conditions, respectively Within each 
set, we select parameters randomly according to the dis- 
tributions described above; see the Appendix for more 
details. 

The model database described above is statistically 
very similar to the database PS-GRB previously used 
to estimate NS-NS and BH-NS merger rates, with one 
critical exception: the simulations used here are ex- 
tensive enough to have a statistically signiflcant num- 
ber of merging BH-BH binaries. Merging BH-BH bi- 
naries are a small fraction of the total BH-BH popu- 
lation and an exceedingly rare consequence of isolated 
binary evolution. For example, the number (n) of BH- 
BH binaries in each of our simulations was always by 
construction 3000 and could be as high as few x 10^ 
(see Figure lll|) . On the contrary, the number (m) 
of potentially merging BH-BH binaries whose delay be- 
tween formation and merger is less than 13.5 Gyr is typi- 
cally 30, multipliedbyaf actorof orderunity; despite our 
efforts, a spiral-galaxy simulation produced only one 
merging BH-BH binary (Figure llip . Further, among 
those few simulated binaries the most massive few can 
be seen to the greatest distance; thus, only a fraction of 
these m binaries significantly impact the average chirp 

mass <^A^c^^/^) [Eq. ®] calculated therefore our es- 
timate of the detection rate. Given the small number 
statistics upon which we build our BH-BH detection rate 
predictions, we very carefully quantify and propagate un- 
certainties into our estimated BH-BH (and also BH-NS, 
NS-NS) detection rates. 

4.1. Model uncertainty and recent literature 
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Our simulations do not explore some known model 
uncertainties that could significantly increase the pre- 
dicted merger rate. As noted previously, our star forming 
condit io ns involve only near-solar metallicity iBulik et alJ 
(pool) : IBelczvnski etall (|2008al [200l . Additionally, 
we assumed that the common-envelope evolution of a 
Hertzsprung gap donor led to stellar m erger, not a com- 
act bi nary. However, as discussed in IBelczvnski et al.l 
2007aD . the alternative and still plausible option will 
increase the BH-BH merger rate by x500. Finally, we 
assume binaries form through isolated evolution alone, 
not allowing for any en hancement in young globular 
prot o-clusters; see, e.g.. lO'Shaughnessv et al.l (l2007b[ ) 
and iSadowski et al.l ()2008i r As mentioned previously, 
because these factors should increase the detection rate, 
our results are the predictions of a conservative model 
family. 

Our simulations also do not sclf-consistently explore 
all model uncertainties corresponding that are compa- 
rable or smaller than the statistical simulation errors 
described extensively below. For example, as the Star- 
Track code has evolved over the extended assembly of 
our archive and this paper, our simulations differ some- 
what from oth er contemporary p apers that employ it. 
Notably, unlike IBelczvnski et al.l (|2007aD , we adopt the 
full Bondi accretion rate during common-envelope evo- 
lution; by trapping more mass in the binary, this gen- 
erally leads to larger detection ratcsQ Also, the Star- 
Track code docs not yet fully explore all uncertainties 
in sin gle star evolution , such as uncertainties in stellar 
radi i (llVver et al.lll999D and due to rotational e ffects (see, 
e.g.. lde Mink et al.l (|2009t ) iMaeder fc Mevnetl pOOl and 
references therein). Nor did we explore modifications 
to the recipes used by StarTrack for single and binary 
evolution. For example, StarTrack models orbital de- 
cay in common-envelope evolution with a single factor 
aA, though in principle A differs between and can be 
calculated for different donor star configurations. For a 
supernova, the initial (pre-fallback) compact object rem- 
nant mass is estimated by t a bulated estimates of the core 
mass; sec Bclc zvnski et al.l ()2008bl ). Compact remnants 
(post-fallback) more massive than 2.5Mq are assumed to 
form black holes. 

Finally, all of our results are sensitive to the total 
number of massive progenitor stars and our estimate for 
the high-mass IMF. While our computationally limited 
exploration cannot give an idealized, fully-marginalized 
theoretical prediction, our approach is extremely useful 

^'^ Though some black hole masses are influenced noticably by 
a change in mass accretion, this change does not always lead to 
drama tic changes in BH-BH m erger- and detection rates. For ex- 
ample, iBeTcz^nskrellin l|2007al ) provide a side- by-side comparison 
of the BH-BH merger rate for two different choices of the accre- 
tion (models A,B in their tables 2,3). They see less than a factor 
2 change between these two extreme limits, much less than the 
typical variation due to the parameters we explore. [See also their 
Table 1 on BH-BH formation channels in these two models.] Given 
our limited computational resources, the fact that BH-BH detec- 
tion rates are ubiquitously dominated by the highest-mass BH-BH 
binaries that merge (not the marginal I3H-BH binaries where one 
formed through AIC), and recongizing that other model parame- 
ters (NS birth mass distribution) and physical inputs (metallicity 
distributions; systematic SFR normalization and model error) will 
perturb our predictions by a comparable amount to the factor of 2 
mentioned above, we are comfortable with neglecting this param- 
eter in this survey of BH-BH detection rates. 



for the reverse problem: how this particular concrete, 
conservative model family can be tested against future 
gravitational wave observations. 

4.2. Properties of all CBC binaries 

Though the calculation was applied only to BH-NS and 
NS-NS mergers, PS-GRB described how to calculate two 
of the three essential ingredients needed to calculate the 
CBC detection rate [Eqs. (^11^ ]: (i) the number of CBC 
merger events per unit mass in progenitors (the mass 
efficiency Ac,q); and (ii) the probability that given a CBC 
progenitor, a merger occurs between t and t+dt since the 
binary's birth as two stars (the delay time distribution 
dPc.a/dt ). For example, as in PS-GRB we estimate 
the mass efficiency A for forming a merging binary of 
type A'(=BH-BH,BH-NS,NS-NS) from the number n of 
binary progenitors of K with 



A = 



n fcut 
N (M) 



(20) 



where N is the total number of binaries simulated, from 
which the n progenitors of K were drawn; (M) is the 
average mass of all possible binary progenitors; and fcut 
is a correction factor accounting for the great many very 
low mass binaries (i.e., with primary mass mi < rric = 
4:Mq) not included in our simulations at all. Expres- 
sions for both (M) and fcut in terms of population syn- 
thesis model parame t ers are provided in Eqs. (1-2) of 
lO'Shaughnessv et al.l ()2007aD . We also estimate the de- 
lay time distribution as before, by smoothing (in logt) 
the n simulated delays t between binary formation and 
merger; see particularly the Appendix of PS-GRB. Fi- 
nally, though PS-GRB does not mention masses, we use 
precisely the same logarithmic smoothing technique to 
estimate dP/dA4c from a set of binaries; see our Ap- 
pendix. 

The procedures described above lead to very reliable 
results for NS-NS and BH-NS binaries, because many 
binaries n and even merging binaries m are present in 
each simulation. However, most simulations have only 
a few BH-BH binaries whose delays between birth and 
merger arc less than the age of the universe. With so few 
merging binaries per simulation, statistical uncertainties 
would severely limit our ability to determine the phys- 
ically relevant portion of delay time (dP/dt) and chirp 
mass {p{A4c)dAic) distributions based on those binaries 
alone. However, since gravitational wave decay depends 
sensitively on the post-supernova orbital parameters of 
a newly-born BH-BH binary, the population of merg- 
ing BH-BH binaries should be very similar to a popula- 
tion with marginally wider orbits but often dramatically 
longer decay timescales. For this reason, in this paper 
we will improve our statistics for dP/dt and p{^Ac) by 
including nearly-merging binaries with delay times t be- 
tween birth and merger less than a cutoff T = 10^ Myr, 
as justified by our studies in the Appendix. 

4.3. Varying the Binary Fraction 

In our population synthesis simulations we systemati- 
cally varied almost all parameters that could significantly 
impact the present-day merger rate. One parameter left 
unchanged in past studies, however, was the binary frac- 
tion /(,, defined as the fraction of stellar systems that 
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are binaries. Without loss of generality, our simulations 
assume all stars form in binaries. The merger rates per 
unit mass implied by a population with a lower binary 
fraction fi, can be related to the merger rates we calculate 
using 

7^(^|^) = 7^(^|^ = l):Ml±||) (21) 

where (q) ~ {1712/ mi) is the average mass ratio of our ini- 
tial stellar population and which is varied between popu- 
lation synthesis models. In particular, because the effect 
of this parameter trivially influences our rate predictions, 
we can analytically incorporate the influence of any prior 
assumptions regarding the binary fraction /t,. 

The true initial binary fraction of stellar populations is 
not known; it is b elieved to be between fb.min = 15% and 
fb,max = 100% (|DuauennoY fc MavoijflQQlt) . For this 
paper we will assume could equally likely take on any 
value in this range. For any given population synthesis 
model, the binary fraction and mass ratio distribution 
multiplicatively influence the observed detection rate by 
a factor 

_ /b(l + (g)) 

which is distributed between Xmm and X ^ 1 (at /t, = 1) 
according to 



p(X)dX 



l + (9> 



dX 



(23) 



4.4. Uncertainties and smoothing : Sampling errors 

Roughly speaking, the relatively small numbers of 
merging binary black holes in any simulation (771 < 100) 
limit our ability to predict the overall detection rate 
log Rd more accurately than l/y^lnlO, because of in- 
accuracies in reproducing both the mass distribution 
dP/dAic and the delay time distribution dP/dt. Addi- 
tionally, the relatively small number of simulations of el- 
liptical and spiral galaxies Ne,s — 300 — 400 provides lit- 
tle chance of discovering high-rate models that are more 
rare than roughly 1/A''. To allow for these two errors, 
we must convolve each individual simulation's merger 
rate with kernels reflecting the poisson uncertainty in- 
troduced into each measurement by the limited number 
of samples. For simplicity, we approximate the relative 
error in log A'^ as normally distributed with standard de- 
viation 1 / \/N In 10. Translating this approximation into 
a PDF, instead of employing a poisson posterior distri- 
bution we adopt a simple gaussian kernel to describe the 
relative error in rate given m and N samples: 

i^sim-A'„(logi?D,l/x/^lnlO) (24) 

i^samp=i^o(logi?D,l/\/AlnlO) (25) 



Ko{x, s) -■ 



:exp(-a;2/2s2) 



(26) 



As noted above, a third uncertainty is the binary frac- 
tion, which we incorporate with the following kernel: 



10-' In 10 



l + (9) 



1- A 



[(g) (10- 



1 



1]^ 



(27) 



for log 0.15 < z < 0. These terms capture the most 
significant simulation-related sources of error. We have 



also explored and included several other potential sources 
of error, including (i) calibration uncertainty in the de- 
tector or mismatch in the waveform model (small, typ- 
ically 10%, multiplied by a factor of order unitvf^: (ii) 
uncertainty in the overall star formation history (a fac- 
tor roughly 2) and in the relative proportions of ellip- 
tical and spiral galaxies (roughly 10%); and (iii) ad- 
ditional sampling errors introduced by the sensitivity 
of the chirp mass average to a few high-mass binaries 
(0.1/ y/rieff / In 10, multiplied by a factor of order unity). 

4.5. Model prediction for CBC rate densities 

To review, using the StarTrack population synthesis 
code, we have explored a representative sample of stellar 
evolution produced by elliptical galaxies {Ne = 206) and 
spiral galaxies (A^5 = 282)0 From these simulations, 
using the tools indicated in § [4] we carefully extracted 
their likely properties (e.g., dP/dt), which in turn we 
convolved with the star formation history of each of the 
two major components of the universe (§ [3]) to gener- 
ate a preferred present-day merger rate per unit volume 
TZk{0) for each model fc, for both elliptical and spiral 
star- forming conditions. 

Though the set of model universe merger rates i?(e,s)fe 
encompasses many of the most significant modeling un- 
certainties, to be more fully comprehensive and to ar- 
rive at a smooth PDF we convolve of the "model uni- 
verse" merger rates rates Rg.p and i?e,/ for p = 1, . . . Ns 
and I = 1 . . . Ne with the uncertainty kernels and bi- 
nary fraction kernel described previously assuming equal 
prior likelihood for each model (i.e., P{Ei) = 1/Ne and 
PiSp) - 1/A^5) 



K = (ATsamp * ATsiin * Kb) 

rip 
Re.l 



(28) 



Pp{\ogRp) = Y,P{Ei)K{\og^) (29) 

(30) 



1=1 

Ns 



p,,(l0gi?,) = ^P(^p)X(l0g-^) 

p—1 

p(l0gi?)-(PeffiPs) (31) 



Figure [5] shows all of the kernels K. In Figure [3] the 
thin {p{\ogR)) and dashed (pe(logi?)) curves correspond 
precisely to the output of this procedure. The bottom 
right panel also shows our estimate for the total merger 

Calibration of the amplitude and phase relationships 
between different measured frequeneies can significantly per- 
turb results. The target calibration amplitude error usu- 
all y cited is 10%, multiplied by a factor of order u nity; 
seelM. Landry (for the LIGO Sci entific Collaboration 11 l|2005l ') and 
lAdhikari et al., (,2003 ) for a detailed survey. 

^■^ We do not calculate BH-NS or NS-NS merger rates for all 
206 elliptical or 282 spiral-galaxy simulations; see Table [T] (avail- 
able in the online version of the Journal) . FoUowring the discussion 
in the Appendix and PS-GRB, we only evaluate BH-NS and NS- 
NS merger rates when they are adequately resolved (n large) and 
(to correct for the bias this introduces) adequately unbiased {nN 
large). Roughly 10% of simulations are not used when evaluat- 
ing the NS-NS or BH-NS merger rate distribution. These models 
form merging BH-BH binaries much more efficiently than BH-NS 
or (occasionally) NS-NS binaries. As discussed in the Appendix, we 
expect observational constraints inevitably rule out the few spiral 
galaxy models involved. 
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Fig. 2. — All of the kernels K versus log X used in this paper, for 
NS-NS (top panel), BH-NS (center panel), and BH-BH (bottom 
panel). As extremely many liinaries are available in each simula- 
tion, the NS-NS simulations' uncertainties are dominated by the 
uncertainty in the binary fraction /j. Because few BH-BH and 
occasionally BH-NS binaries are available in each simulation, how- 
ever, the statistical uncertainty becomes significant. 



rate distribution: a three- fold convolution 

p(log R) = pbh-bh ® Pbh-ns ® Pns-ns (32) 

In the above we assume all simulated population syn- 
thesis models are equally plausible a priori. However, ob- 
servations of single and binary pulsars significantly affect 
our perspective regarding the relative likelihood of dif- 
ferent models. NS-NS merger rates derived from the ob- 
served sample are in fact typically higher than what most 
population models predict for the Milky Way. We adopt 
standard Bayesian methods to determine the posterior 
probability P{Sp\C) for a given spiral-galaxy model Sp 
ip=l...Ns) 

P{Sp\C) cx P{C\Sp)P{Sp) = P{C\Sp)/Ns (33) 

given the observational constraint C on the present-day 
spiral-galaxy contribution to the NS-NS merger rate im- 
plied by present-day observations of pulsars. We adopt 
as a constraint the requirement developed in PS-GRB 
from existing NS-pulsar observations, that the spiral- 
galaxy merger rate density of double neutron stars lie be- 
tween 0.15 and 5.8 Myr"^ Mpc"^ (90% confidence). The 
probabihty P{C\Sp) that the pth spiral-galaxy model 
predicts a rate consistent with observations of dou- 
ble neutron stars is found simply by integrating the 
PDF P{\ogR\Sp) for the rate over the constraint inter- 
vaI3 C = [0.15, 5.8] Myr^^Mpc"^. After renormaliz- 
ing the probabilities P{C\Sp) into P{Sp\C) by requir- 
ing J2p -Pi^pl^) ^ 1' arrive at revised predictions 
for CBC merger rates for BH-BH, BH-NS, and NS-NS 
binaries as in Eqs. (|29ll3ip but with conditional proba- 
bilities: 

Ns 

Ps (log Rs\C) = Y^ P{Sp\C)K {log ) (34) 

P—1 

Convolving the constrained spiral rate density distri- 
bution ps with the unconstrained elliptical density Pe 
leads to our best estimate for the present-day merger 
rate per unit volume, shown as the thick solid curve in 
Figure |31 Our best results favor merger rates between 
2 X 10"''^ - 0.04Mpc^^Myr"^ for BH-BH mergers (90% 

^•^ Strictly, the merger rate distribution function should be cal- 
culated by (i) finding the total merger rate predicted by each model 
(e.g.. Re = Rfybh + Rbhns + Rnsns); (ii) Centering an error kernel 
K on that rate, with n set by the total number of mergers of 
all types; and then (iii) c onv olving the elliptical and spiral distri- 
butions together as Eq. 1311 The approach adopted in the text 
assumes the joint BH-BH, BH-NS, and NS-NS rate distribution is 
uncorrelated, even though weak correlations arc apparent in Table 
[T] across our broad range of parameters (the table is ava ilable in 
the online version of the Journal). As described in ^4.61 this as- 
sumption slightly biases us towards higher rates. On the contrary, 
the three-fold convolution also implicitly introduces independent 
binary fractions for BH-BH, BH-NS, and NS-NS binaries, biasing 
towards lower rates. Empirically, as can be confirmed using Table 
[T] this estimate for p(log R) is in reasonable agreement with the 
summed calculation described above. Finally, since LIGO de- 
tection rates are ubiqutously dominated by the few most massive 
objects that merge — i.e., on Pbh—bh ^ the total detection rate 
does not depend sensitively on this order of evaluation issue. 

Strictly, we should calculate the Bayesian probability with 
some Pl\ogR\0) that describes the merger rate in spirals implied 
by PSR-NS binaries in the Milky Way and uncertainties in the 
star formation rate density of Milky Way galaxies. Our choice of 
a confidence interval based on broadening the the 90% confidence 
limit for PSR-NS binaries simplifies the discussion. 
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confidence), 0.01 - 0.28 Mpc"^ Myr"^ for BH-NS merg- 
ers, and 0.11 - 1.7Mpc~^Myr~^ for NS-NS mergers. 

Our previous papers have similarly constrained binary 
evolution in the Milky Way (PSC2) and in the popu- 
lation of spiral galaxies (PS-GRB). Both interpret ob- 
servations as much tighter constraints than we do here, 
because these papers did not allow for error (e.g., in 
each simulation's prediction and, for the predictions per 
unit volue, in the star formation history of the universe). 
Our bayesian constraints are significantly less restrictive. 
Nonetheless, as seen in Figure |4] the top spiral galaxy 
models (90%) have a kick velocity distribution roughly 
consistent with pulsar observations, similar to that re- 
covered in PSC2 (their Figure 5). Our predictions for 
the posterior NS-NS merger rate are conservative for an- 
other, technical reason: we implement the binary fraction 
as an uncertainty rather than a parameter. We therefore 
cannot rule out the lowest low binary fractions we allow 
ifb = 0.15) in other spiral galaxies, even though a com- 
parison between binary pulsars in the Milky Way and 
our model database disfavors such low values. 

4.6. Discussion 

Two- component universe rarely gives low rates: The 
probability we assign to very low rates depends sensi- 
tively on our implicit assumption that elliptical and spi- 
ral galaxy star forming conditions are uncorrelated. The 
convolution © ps [Eq. (fT8|) ] assigns little probability 
below the sum of the median elliptical and median spi- 
ral merger rate; see, e.g.. Figure El On the contrary, 
if all galaxies have the same undetermined star forming 
conditions, then the relevant merger rate distribution 

p{\ogR) ^^P(g,)^(log f ) (35) 

assigns more proability to the lowest rates; see, e.g., the 
bottom panel of Figure [SI Adopting this implicit prior - 
identical star forming conditions in elliptical and spiral 
galaxies - leads to minimum plausible total BH-BH, BH- 
NS, and NS-NS merger rates 7^.5% roughly 2.4, 7, and 
3.3 times larger than max(72.e_5%, T^.^ 5%), where Tle,5% 
is the 5% probability elliptical-galaxy merger rate and 
similarly; compare, e.g., the top panel of Figure |6] with 
Figure [S] 

Ellipticals have different sta r formation?: As in 
Ide Freitas Pacheco et all (|2006[) . our fiducial calculations 
suggest that elliptical galaxies gave birth to a signifi- 
cant fraction of all present-day compact object mergers, 
particularly BH-BH mergers [Figure [3] . Our results in- 
evitably depend on the amount of and conditions for star 
formation in elliptical versus spiral galaxies. On the one 
hand, we assume most ancient star formation has oc- 
curred in ellipticals. On the other hand, the Salpeter-like 
IMFs {p ~ —2.3) adopted for elliptical galaxies produce 
compact object binaries much more efficiently than the 
much steeper high mass slope {p = —2.7) adopted for 
spiral galaxiesI3 To quantify the impact of each factor 

The steep IMF fo r spiral-like galaxies is motivated by 
IKroupa fc Weianeil 120031 ). which theoretically examined the effec- 
tive single-star IMF produced from continuous formation of open 
clusters. For elliptical galaxies, lacking any definitive suggestion 
that a similar throttled star formation process occurs, we adopt a 
standard Salpeter IMF. 



seperately, in Figure [5] we compare our fiducial BH-BH 
and NS-NS merger rate distributions (thick solid curves) 
with predictions based on only spiral or elliptical star 
forming conditions. Specifically, each new distribution 
still follows from Eq. (|15p [to get TZ — TZg from our 
set of spiral-galaxy simulations] and Eq. (p0| [to cre- 
ate a PDF, based on TZ, an error kernel, and proba- 
bility P{S)]. These comparisons allow us to determine 
the relative impact that three key factors have on our 
predicted results: (i) significantly increasing ancient star 
formation; (ii) requiring ancient star formation produce 
high mass stars more efficiently than present-day star 
formation; and (iii) adopting independent binary evolu- 
tion models for elliptical and spiral galaxies. First, both 
in distribution and on a model-by-model basis, the total 
compact object merger rate due to nearly-steady star 
formation (blue curves, which adopt the "spiral" star 
formation rate p = ps) is less than the amount due to 
all past star formation history (red curves, which adopt 
p ~ Pe + Ps) only by a simulation-dependent factor be- 
tween 1.15-1.3 for NS-NS and between 1.5-2 for BH-BH 
binaries(3 Second, on a modcl-by-modeO and distri- 
bution basis, provided the same amount of star forma- 
tion our elliptical galaxy models produce merging bina- 
ries roughly 1 (NS-NS) or 1 - 10 (BH-BH) times more 
frequently than spiral-only star forming conditions (both 
multiplied by factors of order unity) , entirely due to their 
shallower IMF and broader metallicity distribution; com- 
pare the red and green curves in Figure [S] Third, though 
ellipticals produce merging binaries more efficiently per 
unit mass, fairly few BH-NS and NS-NS binaries remain 
to merge after such a long delay [dP/dt oc 1/t). For BH- 
NS and NS-NS binaries, elliptical galaxies' present-day 
merger rate balances two positive factors - very efficient 
compact binary formation (e.g., ^ 3 times higher, based 
on the IMFs) combined with a high rate of early star 
formation 10 times higher) - against an ^ 30 lower 
merger rate dP/dt via a few Gyr vs ~ lOOAIyr typical 
delay; see the discussion in PS-GRB. For BH-NS and 
NS-NS binaries, this balance roughly leads to compara- 
ble merger rates in elliptical and spiral galaxies To 
summarize, our fiducial unconstrained predictions differ 
from previous predictions based on steady, spiral-galaxy 
star formation by three factors: (i) the addition of an- 
cient star formation, increasing the overall merger rate 
above steady spiral-only by a factor 1.15 — 1.3 for NS- 
NS and 1.5 — 2 for BH-BH; (ii) the conversion of ancient 
star forming conditions from spiral-like to elliptical, in- 
creasing the total present-day merger rate by a factor of 
order unity for BH-NS and NS-NS binaries; and (iii) the 
assumption that elliptical and spiral star forming condi- 
tions are independent, biasing against low total merger 
rates 

The merger rates per unit volume are higher than an extrap- 
olation of the Milky Way's rate to a volume density by a factor of 
roughly [f dP^_NS-Ns/dtp(t)]/p(0)P^^NS-Ns{< 13Myr). 

Though not used or provided here, we have compared several 
pairs of binary evolution models which differ only in their IMF and 
metallicity. 

In fact, the spiral NS-NS merger rate is significantly higher 
than the elliptical NS-NS rate, due to the existence of a population 
of ultracompact binaries not well described by this argument. 

Additionally, we apply observational constraints only to spiral 
galaxies, insuring the retension of elliptical galaxy models with the 
highest BH-BH formation rates. As discussed in ^5.21 the spiral- 
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Fig. 3. — Cumulative probability P{< logi?) of various merger rates per unit volume (units Mpc~^Myr~^), without (thin) and with 
(thick) requiring consistency with the binary pulsar population in the Milky Way. These distributions incorporate uncertainty in the initial 
binary fraction, the unknown population synthesis model parameters, and statistical errors implied by the limited size of our simulations. 
For comparison, the dashed lines show the contribution from ellipticals only; the dotted lines show the contribution from spiral galaxies 
without requiring consistency with the Milky Way. 



Ellipticals "dominate" BH-BH rate: Adopting the same 
binary evolution model but different star formation his- 
tories, elliptical galaxies usually produce more merging 
BH-BH binaries at present than spiral galaxies. In this 
sense, the elliptical galaxy BH-BH merger and detec- 
tion rates "dominate" the total merger rate. However, 
some elliptical galaxy models do yield lower merger rates 
than some spiral galaxy models; neither uhiqutously dom- 
inates. For example. Figure [S] suggests that low but plau- 
sible total BH-NS merger rates require more mergers in 
spirals than ellipticals; high but plausible merger rates 
require the opposite. 

Binary evolution priors: Our predictions depend on the 
relative likelihood of different binary evolution model pa- 
rameters. To facilitate comparison with PSC2 and to 
avoid omitting any vital region of parameter space, we 
allow and treat as equally plausible a wide parameter 
range, including extremely large (and empirically im- 
plausible) supernova kick magnitudes Vkick — 1000 km/s. 
The reader can reproduce our results or explore alter- 
native assumptions by reweighting the data provided in 
Table [1] available in the online version of the Journal. 

galaxy models with frequent-enough NS-NS mergers generally have 
a\ not too small. Very inefficient common envelope evolution is 
required to produce the highest BH-BH merger rates. 



Binary evolution parameters: Finally, as outlined in 
§ 14.11 we fix a number of assumptions which can dra- 
matically influence merger and even detection rates. For 
example, the maximum NS mass m+^rs completely de- 
termines the nature of final compact binaries, as it sub- 
divides a roughly jn+Arg-independent mass distribution 
into three regions: NS-NS, BH-NS, and BH-BH. Evi- 
dently the individual component merger rates depend 
sensitively on m^Ns- However, as relevant gravitational 
wave interferometers cannot easily distinguish between 
two compact objects of comparable mass IMq < m < 
3M0, the overall detection rate described later depends 
only weakly on m-^.NS■ Similarly, while the neutron star 
birth mass distribution can influence the relative pro- 
portions of each type, as well as the observable pulsar 
binary mass distribution, our results should not depend 
significantly on the amount distribution of masses. 

Several other binary evolution parameters not explored 
here have previously been shown to influence merger 
rates significantly, albeit much less significantly than 
the parameters we explored here. First, as noted pre- 
viously in we adopt the full Bondi-Hoyle accretion 
rate onto black holes in common envelope, rather than 
the smaller rate s suggested by r e cent h ydrodynamical 
simulations; cf. IBelczvnski et al.l (|2007al lb[) . A smaller 
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Fig. 4. — Constrained population synthesis parameters for spi- 
ral galaxies: compare to the top right and bottom left panels of 
Figure 5 of PSC2. In each panel, dark points arc included when 
the corresponding model's posterior probability P(Sp|C) is among 
the top 90%; gray points are shown for all other 282 spiral galaxy 
simulations. 



mass accretion rate, if adopted, would reduce the inci- 
dence of accretion induced collapse, marginally chang- 
ing the proportions of BH-BH, BH-NS, and NS-NS bi- 
naries; marginally lower the masses of some merging bi- 
naries, increasing the factor by which common envelope 
evolution tightens orbits; and generally change merger 
and detection rates by a small factor compared with 
the overall range considered [Footnote [10] . Second, we 
examined only isotropic kick distributions. Several au- 
thors have suggested spin-kick alignm ent on theoretica l 
and empirical grounds; see, e.g., Kaloger a et al.l (|2008[ ). 
iKuranov et alj (p009D IWang etal., (,2007,1 and references 
therein. As the same mass transfer that brings binaries 
close enough to merge also aligns their spins and orbit 
before the second supernova, polar kicks are particularly 
efficie nt at disrupting binaries; see iPostnov fc Kuranovl 
(|2008f l. Though polar kicks would dramatically trans- 



form our predictions, the observed PSR-NS merger rate 
in the Milky Way not only corresponds to a high^ not 
low rat^^ but also implies model supernova kick magni- 
tudes comparable to observed pulsar proper motion ve- 
locities [PSC2, as well as our Figure H]. Third, we adopt 
a specific model for mass loss in massive stars; by chang- 
ing the evolution of massive progenitors, alternate stel- 
lar wind models can dramatically modify merger rates 
(cf. iBelczvnski et al.l ()2008af )). Fourth, we adopt a sm- 

"^^ While the range of predicted Milky Way NS-NS merger rates 
depends on the adopted IMF, even a standard Salpeter IMF leads 
to a range of predictions that intersects PSR-NS observations only 
for high rates. 




logR 
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Fig. 5. — Relative likelihood of present-day merger rates R for 
NS-NS binaries (top panel) and BH-BH binaries (bottom panel) 
given five assumptions: (1) only nearly steady spiral galaxy star 
formation produces merging binaries (blue); (2) all galaxies form 
stars like spiral galaxies (red); (3) all galaxies form stars like el- 
liptical galaxies (green); (4) like (l), but using only models that 
reproduce the present-day Milky Way merger rate (blue, dotted); 
(5) like (2), similarly (red, dotted). For comparison, the thick 
black curves shows the fiducial total BH-BH and NS-NS merger 
rate (solid curve includes all spiral models; dotted curve only mod- 
els that reproduce the present-day Milky Way). 



gle, time-independent metallicity for each type of star 
forming region. Binary merger rates can depend sen- 
sitively on metallicity, p articularly through met allicity- 
dcpendent stellar winds (jBelczvnski et al.ll2008al) . A fu- 
ture publication will explore the latter two issues in more 
detail. 

Advanced detectors and redshift- dependent mass evolu- 
tion: For simplicity we have adopted a single time- 
independent mass distributions for each type of merging 
binary: BH-BH, BH-NS, and NS-NS. Looking backward 
to higher redshift, however, the relative proportions of 
each binary type change, as each possesses a different 
delay time distribution dP/dt. In other words, despite 
using "time-independent" mass distributions, the instan- 
taneous total mass distribution, found by gluing these 
three distributions together in proportion to their instan- 
taneous merger rates TZ, will vary with time. Networks 
of advanced detectors and particularly third generation 
detectors will probe a time-dependent mass distribution. 
The authors and collaborators will also address this issue 
in a future paper. 

5. CBC DETECTION RATES AND DETECTION 
PROBABILITY 

5.1. CBC Detection Rates for LICO 



13 




each model (i.e., P{Ei) = 1/Ne and P{Sp) = 1/Ns) 



Fig. 6. — Convolutions illustrated. Top panel; the BH-NS merger 
rates Ti due to ellipticals (dashed), spirals (dotted), and overall 
(solid) assuming uncorrelated elliptical and spiral star forming con- 
ditions, as in the text and Figure [3] The minimum likely merger 
rate overall is significantly larger than the minimum likely due to 
ellipticals or spirals alone. Bottom panel: a gaussian distribution 
p(log a;)cZ log a; (dotted); pffip (solid), analagous to assuming a two- 
component universe where elliptical and spiral galaxy rates are 
independently drawn from p; and the distribution for y = x + x 
(dashed), analagous to a two-component universe with identical 
elliptical and spiral star forming conditions, both drawn from p. 



Using the StarTrack population synthesis code, we 
have explored a representative sample of stellar evolution 
produced by elliptical galaxies {Ne = 206) and spiral 
galaxies {Ns ~ 282). From these simulations, using the 
tools indicated in § |4] we carefully extracted their likely 
properties (e.g., dP/dt), which in turn we convolved with 
the star formation history of each of the two major com- 
ponents of the universe (§[3]). Using our understanding 
of the LIGO range (§ [2]) we can convert these merger 
rate histories into expected single-interferometer LIGO 
detection rates. Finally, though the set of model universe 
detection rates RD,(e.s)k encompasses many of the most 
significant modeling uncertainties, to be more fully com- 
prehensive and to arrive at a smooth PDF we plot in Fig- 
ure [7] the convolution of the "model universe" detection 
rates Rd,s,p and RD,e,i for p = 1, ... iVp and 1 = 1... Ne 
with the uncertainty kernels and binary fraction kernel 
described previously assuming equal prior likelihood for 



K = (i^samp * Ksim * Kb 



Pe{\0gRD,e)=J2P{El)K{\0g 



1 = 1 

Ns 



p,{\ogRn.s)=J2P{Sp)K{\og 

p=i 

p{\0g Rd ) = {Pe ® Ps ) 



RD,e ^ 
RD,e,l 

RD,e N 
RD,e,p 



(36) 
(37) 

(38) 
(39) 



and subsequently combine the elliptical and spiral detec- 
tion rate distributions [Eq. (|18p ]. [Most of these convolu- 
tions can be performed trivially by adding the standard 
deviations of gaussians in quadrature; the remaining two 
convolutions are performed numerically.] As previously 
we estimate the total detection rate via a three-fold con- 
volution [Eq. 132] E3 Though the same basic techniques 
applied above can be and have been applied to BH-NS 
and NS-NS binaries previously (PS-GRB), we now in- 
corporate error propagation, a simulation-by-simulation 
estimate of the relevant chirp mass, and a variable bi- 
nary fraction; the merger rate distributions they provide 
therefore are not proportional to our predictions. 

Figure [7] shows the distribution of CBC detection rates 
predicted by Eq. (fTO)) for individual initial LIGO inter- 
ferometers as part of a network (thin curves) as well as 
the contribution from elliptical galaxies alone (dashed 
curves) As indicated by the close proximity of the el- 
liptical and total merger rate distributions, most BH-BH 
mergers that LIGO detects should be produced in ellip- 
tical galaxies. This elliptical bias arises because BH-BH 
binaries have long delays between their formation and 
eventual merger; their present-day merger rate is more 
easily influenced by ancient star formation. On the other 
hand, the difference between the dashed and thin curves 
in the bottom left panel indicate that the NS-NS merger 
rate must be spiraZ-dominated. 

Because the initial LIGO detectors essentially probe 
different-sized volumes of the local universe for each 
merger type, their detection rate distributions are iden- 
tical to those of any similar-scale single- or multiple- 
IFO network, mod a constant horizontal offset deter- 
mined by the relative increase in volume NS-NS bina- 
ries can be seen [log VArs_Ar5/47r(14Mpc)'^/3]. Similarly, 
because a single advanced LIGO detector's reach to NS- 
NS and BH-NS binaries is often not cosmologically sig- 
nificant, weighting over detector and source orientation, 
our predictions for NS-NS and BH-NS detection rates 
are identical to those for initial LIGO, mod an offset 
[= log(197/14)'^ = 3.44]. However, the typical range of 
a single detector to BH-BH binaries does extend past 
z ~ 0.25, where the star formation rate is noticably 
higher and where cosmological volume and redshift fac- 
tors become significant. But because of the long delay be- 
tween BH-BH merger and progenitor birth, the BH-BH 
merger rate R{t) generally increases much more slowly 

Per the discussion of Footnote 1131 this convolution only ap- 
proximates the true detection rate. In practice, one factor usually 
dominates the detection rate in each component (e.g., for ellipti- 
cals, BH-BH mergers). 

The raw and postprocessed data used to create these figures 
is available from the first author on request. 
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Fig. 7. — Cumulative probability P(< log /Jo) of various event rates for initial (blue; assuming Djj = 31 Mpc) and advanced (green; 
assuming fjj = 445 Mpc) single LIGO interferometers as part of a network, requiring consistency with the binary pulsar population in the 
Milky Way. These distributions incorporate uncertainty in the initial binary fraction, as discussed in the text, as well as uncertainty in the 
many population synthesis model parameters. For comparison the dashed curves show the contribution from ellipticals only. 



with rcdshift than the star formation history. As a re- 
suh, the single-IFO BH-BH detection rate distribution 
is also comparable to but slightly less than what we ex- 
pect after scaling up initial-LIGO results to larger range: 
a detailed calculation suggests a merger detection rate 
between 3 — 300 yr~^ (90% confidence), in good agree- 
ment with the 20 - 500yr-i(£)bns/197Mpc)3 (90% con- 
fidence) expected ft-om cubic rescaling of the initial LIGO 
result For networks of advanced detectors, however, 
cubic scaling breaks down severely; the range is signifi- 
cant enough to require a full cosmological treatment and 
network-sensitivity-beampattcrn averaging. 

5.2. Comparisons with pulsar observations of Galactic 
NS-NS binaries 

Following the previously-discussed procedure to re- 
evaluate the likelihood of our spiral-galaxy population 
synthesis models, we arrive at a distribution of initial 
LIGO detection rates as indicated in Figure [T] Compar- 

■^^ In our simulations, the black hole-black hole rates generally 
are constant or decrease with redshift, due to the long delay be- 
tween binary birth and merger. Also, the "redshifting factor" l + z 
in the detection rate integral implies a slightly smaller comoving 
4- volume swept out by the detector's past light cone per unit time 
at the maximum redshift: J dzdVc/dz/{l + z) ~ Vc/(1 + z). Both 
these factors suggest the BH-BH detection rate at high redshift will 
be lower than a simple cubic extrapolation of the local-universe re- 
sult. 



ing these distributions with the unconstrained results of 
the previous calculation, we arrive at the following con- 
clusions: 

Higher NS-NS detection rate: As seen in Figure |3l in- 
ferences drawn from the double pulsar population about 
the present-day NS-NS merger rate lie on the high end 
of what our models can presently reproduce. By requir- 
ing our models reproduce those observations, our con- 
strained distributions naturally favor higher NS-NS de- 
tection rates. 

Lower BH-BH detection rate from spirals: Compared to 
our previous results (Fig. [7]), the highest binary black 
hole merger and thus detection rates are ruled out in 
spiral galaxies. Physically, the observational constraint 
supports those population synthesis models where one 
particular parameter (aA, related to common-envelope 
evolution) is not exceptionally small {aX > 0.1); see for 
comparison the top left panel of Figure 5 in PSC2. In 
turn, this parameter aX correlates strongly with the BH- 
BH merger rate, with low values of aA being required to 
produce the largest BH-BH merger rates. Thus, because 
observations of NS-NS binaries in the Milky Way favor 
high NS-NS merger rates, they disfavor the highest rates 
of BH-BH mergers in spiral galaxies. 
Observations of the NS-NS merger rate in spiral galaxies 
do not constrain ellipticals, BH-NS well: With the excep- 
tion of NS-NS mergers. Figure [3] demonstrates that most 
constrained predictions for merger rates (and, though not 
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shown here, LIGO's detection rate) agree strikingly well 
with our a priori expectations. These consistent pre- 
dictions should be expected, since observations of spiral 
galaxies cannot constrain dilFerences (e.g., due to mctal- 
licity and IMF) in the binary evolution in elliptical galax- 
ies. In particular, the single most likely source of a gravi- 
tational wave detection in the near future, BH-BH merg- 
ers, are due to their long characteristic ages expected 
to occur primarily in in elliptical galaxies; as a result, 
the present-day BH-BH and total CBC detection rate 
depends only weakly on our ability to rule out certain 
population synthesis parameters in spiral galaxies. Fur- 
thermore, even in spiral galaxies alone, as seen in PSC2 
and references therein, the NS-NS and other CBC merger 
rates are not tightly correlated. Thus despite the fairly 
strong and consistent constraints that double pulsar ob- 
servations imply for binary evolution parameters, as de- 
scribed in detail in PSC2, our overall expectations for 
LIGO's detection rate remain largely unchanged. 

5.3. Net probability of CBC detection 

Factoring in poisson statistics for each model, the net 
prior probability that a single LIGO interferometer 
with a singlc-IFO range to NS-NS inspiral of ZJbns oper- 
ating for time T can make at least one detection is 



-^detect 



-RdT 



)(40) 



where i?bns implicitly enters as a scale factor for Rd 
(oc Dv^). More generally, for a multi-interferometer con- 
figuration with a time-varying network sensitivity, the 
appropriate detection probability can be easily calcu- 
lated from our reference probability by replacing T in the 
above expression by / dt{Dhns{t) / Dhnsref)^': the latter 
quantity characterizes the 4-volume to which the search 
is sensitive. This detection probability can be signifi- 
cantly increased by improvements in the LIGO detector 
or run schedule that improve the detection 4-volume VT. 
Assuming a fixed, orientation-averaged range Z^hns for 
any fraction of the whole LIGO network, the total prob- 
ability of detecting one or more events by that fractional 
network can be approximated by 

Pdetect(> l|14Mpc,T)~0.4 + 0.51og(T/8yr) (41a) 
Pdetect (> 1 1 197 Mpc, T) ~ 0.4 + 0.5 log(T/0.01 y(-)llb) 

for T between lyr and 100 yr or 2 x 10^"^ and 0.1 yr. 
respectively. These approximations are shown in Figure 
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Equivalently, our calculations suggest that a 
gravitational-wave network must have a 4-volume 
sensitivity of 1.3 x lO^Mpc'^yr (roughly equivalent to a 
32 Mpc sphere for 1 year) in order to have a 50% chance 
of detecting one or more events, as 

-^detect (> ll^bns, 1 yr) ^ 0.5 + 1.5 logZ?bns/32 Mp(^42) 

for Z^bns between 20 and 70 Mpc. 

5.4. Comparison to earlier studies 

Our fiducial predictions for compact object merger 
rates [Fig. |3] and gravitational wave detection rates [Fig. 
[7| reflect our uncertainties about the star-formation con- 
ditions and history of the universe. For example, re- 
sults consistent with current observational constraints 




T(yr) 

Fig. 8.— Probability Pdotoct(> ll^^bns.r) [Eq. ^] of de- 
tecting at least one CBC merger in time Tyr by a single in- 
terferometer as part of a network with volume-averaged range 
Dbns = 14 Mpc (solid blue), 32 Mpc (solid red) and 197 Mpc (solid 
green) to double-neutron star inspiral, taking into account infor- 
mation about binary pulsars in the Milky Way. For comparison, 
the das hed curves corresponds to the approximations presented in 
Eq. [41] 



TABLE 1 

Parameters of all binary evolution simulations used in 

THIS paper, along WITH PREDICTED MERGER RATES R, 
DETECTION RATES Rd, AND (FOR SPIRAL SIMULATIONS) POSTERIOR 
PROBABILITIES. THIS TABLE IS AVAILABLE IN THE ONLINE VERSION 

OF THE Journal. 



require compact binary formation in elliptical galaxies 
and therefore disallow the worst-case scenario: mergers 
formed only in spirals by our least-efficient spiral-galaxy 
model. Nonetheless, our broad parameter study encom- 
passes many scenarios and predictions. For 100% bi- 
nary fraction and for each type of compact binary (NS- 
NS, BH-NS, and BH-BH), the most optimistic and pes- 
simistic single-component scenarios in our comprehensive 
table of results (Table [TJ available in the online version 
of the Journal)E3 correspond to merger rates extending 
between lO^"' — 1 Myr^^ Mpc"'^, multiplied by a factor 
of order unity. 

Detailed quantitative comparisons of our results to 
others appearing in the literature must account for 
alternate normalization and star formation history 
choices, which often vary significantly from study to 
study. Our constant star-formation spiral-galaxy pre- 
diction^^ in Table [T] are comparable to Table 4 
in the review article by IPostnqv fc: YungelsonI (j2007D . 



to 



Table I in IBelczvnski et all (I2002D 



and 3 in IBelczvnski et al.1 (|2007al ) 



to Tables 
and to Table 



^* For each simulation used in this study, this table provides 
input parameters; the number of each type of binary appearing 
lOGyr after a starburst (n in the appendix); the number of bina- 
ries predicted to merge within the age of the universe (m in the 
appendix); and, after convolving with a suitable star formation 
rate, the present-day merger and gravitational-wave detection rate 
for each binary type. 

Most predictions in the literature present results only for the 
Milky Way, based on an implicit or explicit steady star formation 
rate assumption. These results can be translated using our present- 
day spiral-galaxy star formation rate (JJJl and the density of Milky 
Way -equivalent galaxi e s, by blue luminosity (~ 0.01 Mpc~'^); see, 
e.e:.. IKopparapu et al.l I I2008I ). 



16 



in iHurlev et alj ( 2002f ). Our constant star- formation 
elliptical-galaxjl^ pr edictions can be com pared to Ta- 
bles 4, 6, and 7 o f 'Voss fc Tauris (2003") and to Ta- 
ble 1 in iPortegies Zw art & Yvmgclson (19981). Our re- 
sults assuming all galaxies have spiral-like star form- 
ing conditions [cf. Fig [5] can be compared with 
iBulik fc Belczvhskil ()2003l ). Finally, our constant and 
time-depen dent elliptical galaxy resu l ts are directly com- 
parable to Ide Freitas Pacheco et al.l (j2006f ) . Even after 
allowing for different inputs, a model-by-model compar- 
ison of our results with others available in the literature 
illustrate physical differences between different authors' 
popul ation models and our own [§ 14. 1| . For example, 
unlike iVoss &: TaurisI ()2003[ ). rather than calculate the 
envelope binding energy fraction A, we fix a single fac- 
tor aX for all binaries in common envelope, leading to 
smaller BH-BH merger rat es for comparable parameters ; 
see, e.g., the discussion in lPostnov fc YungelsonI (I2007D . 
As discussed in unlike IPostnov fc Kuranovl (|2008f ) 
we adopt isotropic kicks, avoiding the least favorable kick 
configurations a nd merger rates; s e e their Figure 3. Un- 
like model C in iBelczvnski et al.l (|2007aD . we assumed 
that the common-envelope evolution of a Hertzsprung 
gap donor led to stellar merger, not a compact binary. 
We apply a similar kick to all NS, rather than assign a dif- 
ferent kick based on its prior (binary) evolution history 
(|Pfahl et all l200l or its supernova mechanism (i.e., a 
low ki ck for electron-cap t ure su pernova; sec Linde n et al.l 
(|2009f ). iKalogera et al.l (I2008D and references therein). 



lO'Shaughnessv et al.l (|2007bD ; lBaneriee et al.l(|2010D and 
references therein). A similar process occur s in nu- 
clear star clusters; see IMiUer fc LaubuF^ ([20091) . Alter- 
natively, young dense clusters could undergo runaway 
stellar collisions, eventually forming an "intermediate- 
mass" black hole M e 100 - lO^Mg, each of which 
can c apture stellar-mass compact objects (jMandel et al.l 
[200l or even form binaries and merge themselves 



Additionally, as described in § 14.11 and § 14.61 we make 
several other specific choices consistent with our previous 
work but not universally adopted in the literature (e.g., 
for the Bondi accretion rate in episodes of hypercriti- 
cal common envelope accretion; for maximum NS mass) . 
However, a detailed exploration of all these correlated 
effects would require coordinated simulations performed 
for the purpose of detailed comparisons; such an under- 
taking is beyond the scope of this paper. Instead, we are 
satisfied that our results are broadly consistent with early 
studies of the Milky Way, particularly given the num- 
ber and parameter survey breadth of the model database 
used here compared with earlier studies by other groups. 
For comparison, in the analysis of our present results 
wc do make detailed comparisons to earlier studies from 
our group (ij 14. Gl and again in our concluding discussion), 
to make the origin of any apparent inconsistencies fully 
transparent. 

Despite providing a distribution, our fiducial results 
conservatively do not account for any other channels 
for forming compact objects besides isolated binary evo- 
lution, even though several observational, semianalytic, 
and numerical studies have demonstrated that interact- 
ing stellar environments can efficiently form compact bi- 
naries. For example, in young very massive stellar clus- 
ters, a high-mass black hole subcluster can mass seg- 
regate and rapidl y form and eject m erging black hole 
binaries (see,e.g., IPortegies Zwart fc McMillan (2000); 

These papers adopt a Salpeter-like IMF but steady star- 
formation rate, a combination not provided in our Tabic [l] For 
a steady star formation rate, the merger rate can be estimated 
from two entries in the table: the ratio of the number of merg- 
ing binaries (m) to the star formation time (T), corrected for the 
mass ratio distribution; see the Appendix of lO'Shaughnessv et al.l 
l|2005l 'l. 



(jFregeau et aI1l2006l; iMandel et al.ll2009[ ). Theory aside, 
observations suggest pulsar and X-ra y binaries ar e highl y 
overabundan t in clusters; see, e.g. . iRasio et al.l (|2000t ): 
IPoolev et all (|2003| ) ; iRansom et"an (|2005l ) and references 
therein. Further, several authors have argued observa- 
tions of short GRBs' hosts, luminosities, and redshift 
distribution support a multi-component origin, as cluster 
dynamics lead t o a different me r ger de l ay-time distribu- 
tion P ( > t): seelHopman et al.l (|2006f ). iSalvaterra et al.l 
(|200l . IGuetta fc Stellal (|2009n and references therein. 
Though present-day massive stellar clusters have little 
mass, if a sufficient fraction of all star formation oc- 
curs through clusters and a sufficient fraction occurs in 
massive clusters, similar highly efficient dynamical pro- 
cesses could produc e many merging comp act binaries of 
aU types; see, e.g., iSadowski et alT (pOM ). To summa- 
rize, both theory and several different types of obser- 
vation suggest that all types of dynamically-produced 
compact binary merger exist. These mergers could occur 
as and potentially orders of magnitude more frequently 
than mergers produced from field binaries; the estimates 
shown in Fig [7] therefore underestimate the total merger 
rate, summed over all channels. Nonetheless, we antici- 
pate that future observations might distinguish between 
the relative contributions of each component. Assuming 
even rough separation, our results determine how and 
how precisely future measurements of compact object 
merger rates (through gravitational waves, short GRBs, 
or galactic observations) will constrain binary evolution. 

6. CONCLUSIONS AND DISCUSSION 

In this paper, we provide the first estimate of the CBC 
detection rate for initial and advanced LIGO which ac- 
counts for (i) star formation in both elliptical and spiral 
galaxies as a function of redshift, (ii) binary black holes, 
(iii) a range of plausible binary fractions, and most sig- 
nificantly (iv) a large range of plausible binary evolution 
scenarios. Our results strongly indicate that that exist- 
ing gravitational wave detectors are at the threshold of 
detection; that moderate improvements in these detec- 
tors, such as the enhanced LIGO upgrade of 2009, could 
plausibly lead to detecting at least one binary merger 
in a year of operation; that advanced LIGO detectors 
should be reasonably expected to detect inspirals over 
1-2 yr of operation; and that even the absence of a de- 
tection by advanced interferometers would very signifi- 
cantly constrain the set of model parameters that could 
be consistent with observations. 

Our predicted detection rates are somewhat higher 
than those presented in PSC2 because we include ellipti- 
cal galaxies: mergers in elliptical galaxies likely dom- 
inate our estimates of the BH-BH and BH-NS event 
rates. However, our previous estimates effectively as- 
sumed a high present-day star formation rate per unit 
volume: an assumed Milky Way-equivalent galaxy den- 
sity of 0.01/ Mpc"^ and star formation rate per spiral 
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galaxy of 3.5MQ/yr correspond to a present-day star 
formation rate three times higher than observed. For 
this reason, a naive product of the merger rates per 
Milky Way galaxy drawn from PSC2's Figure 6 leads to 
preferred merger rates per unit volume that agree well 
our preferred values. For example, our current results 
and this naive rescaling of PSC2 results both favor lo- 
cal merger rate densities of 10~^ Mpc~^ Myr~^ for black 
hole-neutron star binaries and 3 x 10~^ Mpc~^ Myr"^ for 
binary neutron stars. Correcting for the star formation 
history normalization, however, our calculations suggest 
three times as many mergers per volume for BH-NS and 
NS-NS than in PSC2. Similarly, our predicted detection 
rates are higher than those implied by Figure 7 in PSC2, 
correcting for the previously mentioned star formation 
rate bias (reducing the rates shown in that figure by 3) 
and rescaling to a single-interferometer rather than net- 
work range (reducing the rate by xlO). 

Limiting attention to elliptical- plus spiral-galaxy mod- 
els that include only those spiral galaxy history mod- 
els that reproduce current observations of double pul- 
sars, we find that a single detector in the initial 
LIGO network should be sensitive to 2.4 x 10~^— 0.46 
(Pc/8)"^(-Db„s/14Mpc)^ mergers per year (90% confi- 
dence level), where we sum over all merger types to 
produce an overall detection rate, where = 8 is the 
threshold SNR of the search and Di,ns is the radius 
of the effective volume inside of which a single detec- 
tor in the appropriate array could observe the inspi- 
ral of two 1.4M0 neutron stars. We estimate that the 
LIGO detector network has a probability Pdetcct(> 1) > 
0.4 -I- 0.5 log(y/V'c)(T/8 yr) of detecting a merger duing 
each year of operation, where V = 47rZ?bns'^/3 is the 
volume to which a multidetector search for binary neu- 
tron star inspiral is sensitive, Vc is the volume inside a 
14Mpc radius sphere, and 1/3 < (V^/K)(T/yr) < 8. 
For example, the probability Pdctcct of detecting one 
or more mergers can be approximated by (i) Pdctcct — 
0.4 + 0.51og(T/0.01yr), assuming Z^bns = 197Mpc and 
it operates for T years, for T between 2 days and 0.1 yr); 
or by (ii) Fdctcct ^ 0.5-1-1. Slog Dbns/32Mpc, for one year 
of operation and for -Dbns between 20 and 70 Mpc. 

In this paper we have continued to employ approx- 
imate waveform and detection models to estimate the 
sensitivity of ground-based gravitational wave detectors 
to m oderate- mass (M :^ 1 Mp,) BH- BH mergers (see, 
e.g., Flanagan fc Hughe j (|1998l ) and iBuonanno et al.l 
(|2003| ) for a review)^ Advanced detectors, however, 
will be sensitive at cosmological distances to the later 
phases of high-mass mergers, a complex strong-field and 
strongly-precessing regime which is only beginning to be 
thoroughly understood through numerical simulations. 
Progress in numerical waveform modeling and search pi- 
plcline analysis is required to better understand how of- 
ten these candidates could be identified, especially for 
spinning systems, as expected in nature. 

6.1. Anticipated constraints from future observations 

Despite the existence of two competing scenarios for 
compact binary formation, in many cases the electro- 
magnetic and gravitational wave signatures of bina- 
ries formed dynamically and in isolation can be distin- 
guished. For example, black hole binaries formed dynam- 



ically should have random spins and therefore (gener- 
ically) beating gravitational inspiral waveforms, while 
isolated binary evolution models favor spin-orbit align- 
ment. Once binaries with isolated progenitors are iden- 
tified, their gravitational wave signals will encode both 
the number and properties of binaries formed through 
isolated evolution, information that in turn can tightly 
constrain isolated binary progenitor parameters. To use 
an extreme example, future BH-BH observations could 
imply a gravitational- wave detection rate Rd consistent 
with only the highest event rates shown in Fig. [71 In the 
context of the models presented here, such an observation 
would imply (i) that elliptical galaxies host most BH-BH 
binaries; (ii) that massive stars have a high binary frac- 
tion; (iii) that a certain phenomenological parameter is 
very small {aX < 0.1; H5.2^ : and thus (iv) since Milky 
Way observations of binary pulsars exclude low aX, that 
isolated binaries formed in elliptical galaxies evolve dif- 
ferently to those in spirals. However, by analogy with 
the insights provided from observations of isolated Galac- 
tic pulsars, typical detection rates will be compatible 
with many plausible progenitor models; see PSC2 and 
[0_'Sliauglmessy ct al. (2005). Rather than constraining 
single parameters tightly, the information provided by 
the first few observations will be encoded only in high- 
order correlations among many model parameters; see, 
e.g.. Fig. 5 in PSC2 and the hig hly correlated merger 
rate fi t presented in Appendix B in lO'Shaughnessv et al.l 
(120051 ). Finally, what we can learn about binary evolu- 
tion depends on the merger detections that nature pro- 
vides. Because correlations between parameters change 
significantly across the model space, we cannot describe 
posterior uncertainties once and for all. 

If neutron star mergers occasionally have electro- 
magnetic counterparts (i.e., short GRBs), then well- 
localized electromagnetic counterparts to binary neu- 
tron star merg ers, such as (promp t or orphan ) shor t 
GRB afterglow s [Gehrels et all (l2009l): iRau et al.l (|2009l): 
IBloom fc et ail (120091) or radio bursts, allow optical 
follow-up studies of its host galaxy's star-forming and 
chemical enrichment history. Combined with the po- 
sition of the burst on the galaxy (kick) and, if within 
range of gravitational wave detectors, the binary compo- 
nent masses, this host information can tightly limit iso- 
lated binary progenitor models; see, e.g., PS-GRB. For 
example, these observations could provide unique direct 
probes on the distributi on of de l ays be tween bina ry birth 
and m e rger; see, e.g.. iNakarl (l2007h. IProchask a et al.l 
(|2006[) . iZheng fc Ramirez-Ruid (I2007D . iRhoadsi ()20100 
and references therein. Better still, offsets from host 
galaxies might provide an independent constraint on de- 
lay times (and on supe r nova kick ve l ocities ); see, e.g., 
iBelczvnski eTaD (p006l ). iTroia et al.l (pOOl ) and refer- 
ences therein. Unfortunately, purely electromagnetic ob- 
servations of short GRBs cannot distinguish between 
multiple source populations, such as cluster-formed 
mergers or even non-merger-powered bursts (see, e.g. 
iVirgih et al1l2009t iChapman et al.l 120081 ). Worse, given 
plausible kick offsets, purely electromagne tic observa- 
tions may not isolate the appropriate host (jZemp et al.l 
I2009D . Fortunately, if existing space-based and ground- 
based resources are applied in the advanced LIGO epoch, 
coincident gravitational- wave and electromagnetic obser- 
vations can resolve these ambiguities. 
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Finally, as always with binary evolution, any further 
observations pertaining to binary stars can improve our 
understanding and limit model-building freedom, such as 
SNIa rates; galactic pulsar and X-ray binaries; and even 
ga lactic WD binarie s [see, e.g., the extensive discussion 
in lNelemans fc et~a. I (l200l ]. 



19 



The authors appreciate helpful comments received 
from Ilya Mandel, Alberto Vecchio, Patrick Brady, Chad 
Hanna, Steve Fairhurst, Tania Regimbau, Alan Wein- 
stein, and all the members of LIGO inspiral search 
group. Grants R.O. was supported by National Sci- 
ence Foundation awards PHY 06 -53462 and the Center 
for Gravitational Wave Physics. The Center for Gravita- 



tional Wave Physics is supported by the George A. and 
Margaret M. Downsbrough Endowment and by the Na- 
tional Science Foundation under cooperative agreement 
PHY 01-14375. VK was supported by NSF grant PHY- 
0653321. KB acknowledges the support from the Polish 
Ministry of Science and Higher Education (MSHE) grant 
N N203 302835 (2008-2011). 



REFERENCES 



Abbott et al. (The LIGO Scientific Collaboration). 2003, 

(gr-qc/0308043) [URL] 
— . 2006. Phys. Rev. D, ^3, 062001 

[ADS] 

— . 2008, Phys. Rev. D, 77, 062002 "[URL]' 

Abbott et al. (The LIGO Scientific Collaboration), B. 2009, 
Phys. Rev. D, 79, 122001 [ADS] [ADS] 

Acernese, P., Amico, P., Alshourbagy, ivi., Antonucci, P., Aoudia, 
S., Avino, S., Babusci, D., Ballardin, G., Barone, P., Barsotti, 
L., Barsuglia, M., Beauville, P., Bigotta, S., Birindclli, S., 
Bizouard, M. A., Boccara, G., IBondu, P., Bosi, L., Bradaschia, 

C. Braccini, S., Brillet, A., Brisson, V., Brocco, L., Buskulic, 

D. , Calloni, E., Campagna, E., Cavalier, P., Cavalieri, R., 
Cella, G., Gesarini, E., Chassande-Mottin, E., Corda, C., 
Cottone, P., Clapson, A.-C, Cleva, P., Coulon, J. -P., Cuoco, 

E. , Dari, A., Dattilo, V., Davier, M., DeRosa, R., Di Fiore, L., 
Di Virgilio, A., Dujardin, B., Eleuteri, A., Enard, D., Perrante, 
I., Pidecaro, P., Piori, I., Plaminio, R., Pournier, J.-D., 
Francois, O., Frasca, S., Frasconi, P., Freise, A., Gammaitoni, 
L., Garufi, P., Gennai, A., Giazotto, A., Giordano, G., 
Giordano, L., Gouaty, R., Grosjean, D., Guidi, G., Hebri, S., 
Heitmann, H., Hello, P., HoUoway, L., Karkar, S., 
Kreckelbergh, S., La Penna, P., Laval, M., Leroy, N., Letendre, 
N., Lorenzini, M., Loriette, V., Loupias, M., Losurdo, G., 
Mackowski, J.-M., Majorana, E., Man, C. N., Mantovani, M., 
Marchesoni, F., Marion, P., Marque, J., Martelli, P., Masserot, 
A., Mazzoni, M., Milano, L., Moins, C, Moreau, J., Morgado, 
N., Mours, B., Pai, A., Palomba, G., Paoletti, P., Pardi, S., 
Pasqualetti, A., Passaquieti, R., PassucUo, D., Pcrniola, B., 
Piergiovanni, P., Pinard, L., Poggiani, R., Punturo, M., Puppo, 
P., Qipiani, K., Rapagnani, P., Reita, V., Remillicux, A., Ricci, 
P., Ricciardi, I., Ruggi, P., Russo, G., Solimeno, S., Spallicci, 
A., Stanga, R., Taddei, R., Tonelli, M., Toncelli, A., Tournefier, 
E., Travasso. P., Vajente, G., Verkindt, D., Vetrano, P., Vicero, 
A., Vinet, J.'-Y., Vocca, H., Yvert, M., & Zhang, Z. 2006, 
Classical and Quantum Gravity, 23, 635 

r[ADS] 

Adhikari, R., Gonzalez, G., Landry, M., O'Reilly, B., & the LIGO 
Scientific GoUaboration. 2003, Classical and Quantum Gravity, 
20. 903 

[ADS] 

Ajith, Jr"., Hannam, M., Husa, S., Chen, Y., Bruegmann, B., 
Dorband, N., Mueller, D., Ohme, F., PoUney, D., Rcisswig. C, 
Santamaria, L., & Seller, J. 2009, (arXiv:0909.2867) [URL] 

Anderson, W. G., Whelan, J. T., Brady, P. R., Creighton, J. 
D. E., Chin, D.', & Riles, K. 2001, Beam Pattern Response 
Functions and Times of Arrival for Earthbound Interferometers, 
Tech. Rep. LIGO-TOIOIIO-OO-Z, LIGO Project [URL[| 

Apostolatos, T. A., Cutler. C, Sussman. G. J., & ihorne, K. S. 
1994. Phys. Rev. D. 49, 6274 
[ADS] 

Banenee. S.. Baumgardt, H., & Kroupa, P. 2010, MNRAS, 402, 

371' [URL] [ADS] 
Belczynski, K., euiik, T., Fryer, C. L., Ruiter, A., Vink, J. S., & 

Hurlev, J. R. 2009, ArXiv e-prints 

[ADS] 

Belczynski, K., Hartmann, D. H., Fryer, C. L., Holz, D. E., & 
O'Shea, B. 2008a, ArXiv e-prints 

[ADS] 

Belczvnski, K., Kalogcra, V., & Bulik, T. 2002, ApJ, 572, 407 
[ADS] 

Belczynski, K., Kalogcra, V., Rasio, P., Taam. R.. Zezas. A.. 

Maccaronc, T., & Ivanova, N. 2008b. ApJS,' 174, 223 [URL]] 

[ADS] 

Belczynski, K., Perna, R., Bulik, T., Kalogera, V., Ivanova, N., & 

Lamb. D. Q. 2006, ApJ, 648, 1110 

[ADS] 

Belczynski, K., Taam, R. E., Kalogcra, V., Rasio, P. A., & Bulik, 
T. 2007a, ApJ, 662, 504-|- 
[ADS] 

Belczynski, K., Taam, R. E., Rantsiou, E., & Sluvs. M. v. d. 

2007b, submitted to ApJ (astro-ph/0703131) [tJRL] 
Bloom. J. S. & et al. 2009, ArXiv e-prints ' 

[ADS] 



Brown, D. 2004, PhD. Thesis, U. Wisconsin at Milwaukee I [URLjl 

Brown. G. E. & Bethe, H. A. 1994, ApJ, 423, 659 

'[ADS] 

Buhk. T. & Belczynski, K. 2003, ApJ, 589, L37 
[ADS] 

Buhk, i., Belczynski. K.. & Prestwich, A. 2008, 

(arXiv:0803.3516) [URL] 
Bulik, T., Gondek-R osihsK a, D., & Belczynski, K. 2004, MNRAS, 

352, 1372 

[ADS] 



Buonanno, A., Chen, Y., & Vallisneri, M. 2003, 67, 104025 1 [URL] I 
Buonanno, A., Chen, Y., & Vallisneri, M. 2003, Phys. Rev. D, faV, 
104025 

'[ADS] 

Gnapman, R., Priddev. R. S., & Tanvir, N. R. 2008, 

(arXiv:0802.0008) [URL] 
Cutler. C. & Flanagan, E. 1994, Phys. Rev. D, 49, 2658 

[ADS] 

de treitas Pacheco, J. A., Regimbau, T., Vincent, S., & Spallicci, 
A. 2006. International Journal of Modern Physics D, 15, 235 
'[ADS] [ADS] 

de Mink. s. t^., Cantiello, M., Langer, N., Pols, O. R.. Brott, I., & 
Yoon, S.-C. 2009, A&A, 497, 243 
[ADS] 

Duquennoy, A. & Mayor, M. 1991. A&A, 248, 485 
[ADS] 

Pah. S. M., Chandar, R., & Whitmore. B. C. 2005, ApJ, 631, L133 

[ADS] 

Finn. L. S. & Chernoff, D. P. 1993, Phys. Rev. D, 47, 2198 

[ADS] 

Flanagan, E. E. & Hughes, S. A. 1998, Phys. Rev. D, 57, 4535 

[ADS] 

Pregeau, J. M., Larson, S. L., Miller, M. C., O'Shaughnessy, R., 
& Rasio, P. A. 2006, ApJ, 646, L135 [ADS] [ADS] 

Prver. C, Burrows, A., & Benz, W. 1998, Apj, 496, 333 
[ADS] 

Fryer, (J. L., Woosley, S. E., & Hartmann, D. H. 1999, ApJ, 526, 
152 

[ADS] 

Gallagher. HI, J. S., Hunter, D. A., & Tutukov, A. V. 1984, ApJ, 
284, 544 
[ADS] 

Gehreis, N., Ramirez-Ruiz, E., & Pox. D. B. 2009, ARA&A. 47, 
567 

J [ADS] 

Giietta. D. & SteUa. L. 2009, A&A. 498. 329 

[ADS] 

Hopman, C Guetta, D.. Waxman, E., & Portegies Zwart, S. 
2006. ApJ, 643, L91 
[[ADS] 

Hurlev. J. R., Tout, C. A., & Pols, O. R. 2002, MNRAS, 329, 897 
[ADS] 

Kalogcra, V., Valsecchi, P., & Willems, B. 2008, in American 
Institute of Physics Gonference Series, Vol. 983, 40 Years of 
Pulsars: Millisecond Pulsars, Magnetars and More, ed. 
C. Bassa, Z. Wang, A. Gumming, & V. M. Kaspi, 433-441 
Kennicutt, Robert C, J. 1998. Annual Review of Astronomy and 

Astrophysics, 36, 189 [URXJ 
Kopparapu, R. K.. Hanna, ^. R., Kalogera, V., O'Shaughnessy, 
R.. Gonzalez. G.. Brady. P. R.. & Fairhurst, S. 2008, ApJ, 675, 
1459 [URL] [ADS] 



006, 



Kroupa, P. & Weidner, C. 2003. ApJ. 598. 1076 [URL] I 
Kulczycki, K., Bulik, T., Belczynski, K., & Ruda Tt, B. 2 i 
A&A, 459, 1001 

•[ADS] 

Kuranov, A. G., Popov, S. B., & Postnov, K. A. 2009, MNRAS, 
395. 2087 

[ADS] 

Li, ij., Zittang, F., & Han, Z. 2006, accepted for publication to 

ChJAA (astro-ph/0605610) [URL] 
Linden, T., Sepinsky, J. P., Kalogera, V., & Belczynski, K. 2009, 

AdJ, 699, 1573 

|[ADS]| 



20 



M. Landry (for the LIGO Scientific Collaboration). 2005, 
Classical and Quantum Gravity, 22, S985 [URL] 

Maeder. A. & Meynet, G. 2001, A&A, 373, 005 " 
[ADS] 

Mandel, L, Brown, D. A.. Gair, J. R., & Miller. M. C. 2008, ApJ, 
681. 1431 

[ADS] 

Mandel. I., Gair. J. R., & Miller. M. C. 2009, ArXiv c-prints 
[ADS] 

Miller. M. C. & Lauburg, V. M. 2009, ApJ, 692, 917 

[ADS] 

Nagamme, K.. Os triker. J. P.. Fukugita. M., & Con, R. 2006, 

ApJ, 653, 881 [URL]J[ADS]| 
Nakar. E. 2007, Phys. Rep., 44 2, 166 

[ADS] 

Neiemans, G. & et al. 2009, in Astronomy, Vol. 2010, AGS Stars 

and Related Phenomenastro2010; The Astronomy and 

Astrophysics Decadal Survey, 221 — h 
O'Shaughnessv. R.. Kalogera, V., & Belczynski, C. 2008a, ApJ, 

675, 566+ [URL] 
O'Shaughnessv. rt., Kalogera, V.. & Belczynski, K. 2007a, ApJ, 

667, 1048 [URL] 
O'Shaughnossy, K., Kim, C, Fragos. T.. Kaloeera. V.. & 

Belczynski, K. 2005, ApJ, 633, 1076 |ADS] [ADS] 
O'Shaughnessy, R., Kim. C. K alogerarV.. & tseiczynski, K. 

2008b, ApJ, 672, 479 [UR^] 
O'Shaughnessy, R., Kopparapu, R., & Belczynski, K. 2008c, 

submitted to ApJ; arXiv:0812.0591 [URL] 
O'Shaughnessy, R., O'Leary. R.. & Rasio, . A. 2007b, 

Phys. Rev. D, 76, 061504 [URL] 
Peters. P. 1964, Physical Review. 136. 1224 

[AD'S] 

Plahi, li., Rappaport, S., Podsiadlowski, P., & Spruit, H. 2002, 
AdJ. 574, 364 
[ADS] 

Pooicy, D., Lcwin, W. H. G., Anderson, S. F., Baumgardt, H., 
Filippenko, A. V., Gaensler, B. M., Homer, L., Hut, P., Kaspi, 
V. M., Makino, J., Margon, B., McMillan, S., Portegies Zwart, 
S.. van der Klis, M., & Verbunt, F. 2003, ApJ, 591, L131 
[ADS] 

Portegies Zwart, S. F. & McMillan, S. L. W. 2000, Astrophysical 

Journal. 528, L17 [URL] 
Portegies Zwart, S. t . k Yungclson, L. R. 1998, A&A, 332, 173 

[ADS] 

Postnov, K. & Yungclson, L. 2007, ( astro- ph/0701059) |[URL]| 
Postnov. K. A. & Kuranov, A. G. 2008, MNRAS, 384, 1393—^ 
[URL] [ADS] 



Prochaska, J. X., Bloom, J. S., Chen, H.-W., Foley, R. J., Perley, 
D. A., Ramirez-Ruiz, E., Granot, J., Lee, W. H., Pooley, D., 
Alatalo, K., Hurley, K., Cooper, M. C, Dupree, A. K., Gerke, 

B. F., Hansen, B. M. S., Kalirai, J. S., Newman, J. A., Rich, 
R. M., Richer, H., Stanford, S. A., Stern, D., & van Breugel, 
W. J. M. 2006, ApJ, 642, 989 

II [ADS] 

Ransom, S. M., Hessels, J. W. T., Stairs, L H., Freire, P. C. C, 
Camilo, F., Kaspi, V. M.. & Kaplan, D. L. 2005. Science, 307, 
892 

[ADS] 

Rasio. F. A., Pfahl, E. D., & Rappaport, S. 2000, ApJ, 532, L47 
[ADS] 

Rau. A.. Kulkarni, S. R., Law, N. M., Bloom, J. S., Ciardi, D., 
Djorgovski, G. S., Fox, D. B., Gal- Yam, A., Grillmair, C. C, 
Kasliwal, M. M., Nugent, P. E., Ofek, E. O., Quimby, R. M., 
Reach, W. T., Shara, M.. Bildsten, L.. Cenko, S. B., Drake, 
A. J., Filippenko, A. V., Helfand, D. J., Helou, G., Howell, 
D. A., Poznanski, D., & Sullivan, M. 2009, (arXiv:0906.5355) 
TURL] 

Rnoads. J. E. 2010, ApJ, 709, 664 
[ADS] 

SadowsKi. A.. Belczynski, K., Bulik, T.. Ivanova, N., Rasio, F. A., 
& O'Shaughnessy, R. 2008, ApJ, 676', 1162 
[ADS] 

Salvaterra, R., Cerutti, A., Chincarini, G., C olpi. M .. Guidorzi, 

C, & Romano, P. 2007, (arXiv: 0710.30 99) [URL]| 
The VIRGO project website. ???? |[URL]| 

Troja, E., King, A. R., O'Brien, P. T., Lyons, N., & Cusumano, 
G. 2008, MNRAS, 385, LIO 
f[ADS] 

Virgih, F. J., Zhang, B.. O'Brien, P., & Troja, E. 2009, 

(arXiv:0909.1850) [URL] 
Voss, R. & Tauris, 1 . M. 2003, Monthlv Notices of the Royal 

Astronomical Society 342, 1169 [URL] 
Wang. C. Lai, D., & Han, J. L. 200^ ApJ, 656, 399 

[ADS] 

Zemo. M., Ramirez-Ruiz. E., & Diemand. J. 2009. ApJ, 705, L186 
[ADS] 

Zhang. Q. & Fall. S. M. 1999, Astrophysical Journal, 527, L81 
[URL] 

Zheng, Z. & Ramirez- Ruiz, E. 2007, ApJ, 665, 1220 1 [URL] | [ADS] | 



APPENDIX 



RECONSTRUCTING THE IVIASS AND TIME DISTRIBUTIONS 

Selecting binaries for mass estimation: Wc prefer to smooth dP/dt over a shorter timcscalc than the maximum allowed 
by evolution of isolated BH-BH binaries. On physical grounds, double BH binaries whose delay times are comparable 
to or smaller than the age of the universe are more likely to be similar to merging BH-BH binaries (the focus of our 
study) than wide binaries that have never interacted. To quantify the point at which a the transition occurs and 
therefore the timescale inside of which BH-BH binaries are comparatively similar, we calculate the "average" chirp 
mass of binaries Aic{t) merging within a factor 10 of time t: 



Mc{t)-- 



Si:0.1t 



<tj<10t-'^'cj 



X]j:t0.1<tj<tlO ^ 



6/15 



(Al) 



Figure ini shows the results of this calculation. Based on the delay time distribution Pm{< t) and average chirp mass of 
similar binaries M.c{t)^ a sharp distinction exists between binaries with t < 10^ Myr and wider binaries. These massive 
wide binaries have orbits largely unaffected by binary evolution and have a delay time distribution dP/dt (x 1/t 
largely determined by their initial orbital configuration. On the contrary, tighter and less massive binaries have been 
significantly perturbed by binary evolution. Within each simulation, they appear to be drawn from a similar chirp 
mass distribution (allowing for significant sampling- induced fluctuations in Aic{t)) 

Predictions and smoothing: Physical timescales From the above similarity argument, we suspect the ?^e// BH-BH 
binaries with delay times times < 10'^ Myr are similar to the set of binaries with delay times < 13.5 Gyr. However, 
even with those binaries we often have only ^ 50 with which to estimate p(A^c); the chirp mass distribution must be 
estimated with smoothing. We use the same approach as is used for dP/dt in PS-GRB to estimate the chirp mass 
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Fig. 9. — Top panels: Cumulative probabilities Pm{< t) that a BH-BH binary will merge in time less than t, for ten randomly-chosen 
population synthesis models that satisfy the limitations discussed in the text, given spiral (left) and elliptical (right) star forming conditions. 
A vertical dashed line indicates the age of the universe. Given the enormous sample sizes involved - by construction, each sample has at 
least 3000 binaries - these distributions are expected accurate to within 0.03 almost everywhere (with 90% probability), barring the fine 
details near the short-time and long-time limits. Bottom panels: Average chirp mass Mc(t) [Eq. (|AT}] of binaries merging between O.lt 
and lot, versus time. Usually, the binaries that merge within 100 Gyr after their birth have a chirp mass lower than IOA/q, significantly 
lower than the average chirp mass of more separated holes. 
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Fig. 10. — Detailed view of the cumulative probabilities Pm(< t) for 4 randomly selected spiral (left) and elliptical (right) models versus 
logf/Myr; compare with Figure[9] Solid curves indicate the fit; points indicate the raw cumulative Pm{< t)- The horizontal dash ed lines 
indicate the smallest and largest possibilites for logP logl/n, the smallest resolvable cumulative probability; see Figure im and the 
associated text for a discussion of its range and relation to m (i.e., the largest 1/n occurs in only a handful of simulations, most of which 
have many merging binaries m S> 1 and therefore a wcU-resolvcd P(< t)). Our estimate reproduces Pm(< t) reliably at t ~ lOOMyr, the 
smallest scale on which p(i) varies for the redshift range of interest. Though our fits necessarily fail to reliably recover dPm/dt on ve ry 
short timescales, uncertainty in integrals of dPm/dt over longer timescalcs average out (i.e., in merger rate estimates 7?.(t) using Eq. I15I I 
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Fig. 11. — Left plot: For spiral (top panels) and elliptical (bottom) archives, a scatter plot of the number of BH-BH binaries n seen in a 
simulation of N progenitor binaries (each progenitor component having mass greater than 4Mq ; higher primary masses are drawn from the 
brokcn-Kroupa IMF). The two solid lines show the cutoffs n > 3000 and nN > 1 X 10^ imposed to insure data quality and reduce sampling 
bias. Simulations used in this paper, shown in blue, must lie above and to the right of these cutoffs. Right plot: For the simulations above 
and to the right of the lines shown to left, scatter plot of the total number of bound BH-BH binaries remaining in the simulation (n) versus 
the number m of potentially-merging binaries (i.e., with a delay between birth and merger less than 13.5 Gyr). All selected simulations 
{uN > 1 X 10^) have at least one merger (m > 0). 



distribution: 



p(x J = ; y ' ,-iiosM.-i^,,r/2is,r (A2) 

_ [(maxfc^m.fc) - (mmklm^k)] „x 
1.25^71777 

where the lm,k = logio(-^c/.^0) for fc = 1 . . . n^ff are the logarithms of each binaries' chirp mass for the ?^e// binaries 
with delay times < 10^ Myr. 

Absent correlations: We identify very few significant correlations on the timescales of interest between log A^c a-nd logt. 
For example, the correlation coefficient between logTMcfc and \ogtk for each simulation's fie// close BH-BH binaries 
is roughly consistent with sampling error (cx 1/ ^Jn^JJ). For BH-NS and particularly NS-NS, nominally significant 
correlation coefficients occur. For example, some NS-NS (A^c,i) distributions concentrate a fraction of mergers on 
and exclude a region below a critical t{M.c) curve, associated with mass-transfer- limited ultracompact mergers on 
times t < 100 Myr. Other distributions show no apparent structure, but through large ri^ff and therefore extremely 
narrow mass distributions yield nominally significant correlation. Our manual inspection of all these outliers suggests 
no structure significant enough to impact merger rates (i.e., associated both with significant change in A^c and long 
star formation timescales). On the contrary, as seen in Figure [9l the joint distribution p(A^c,i) varies significantly 
between simulations, as the distribution in each factor [Mc and t] changes significantly. The authors can provide all 
binary properties [(mi, m2, t)] for each simulation on request. 

SAMPLING LIMITATIONS AND ARCHIVE SELECTION 

Both sets of simulations are highly heterogeneous: the number of binaries simulated (N), the number of bound BH- 
BH binaries remaining at the end of stellar evolution (n), and the number of those binaries which merge within 13.5 Gyr 
(m) all span several orders of magnitude, as shown in Figure 1111 Evidently, only for those simulations with many 
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BH-BH binaries {n ^ 1000) - particularly with several merging binaries {m 3> 1) - can we reliably extract relevant 
physics, such as the present-day merger rate and characteristic chirp mass of BH-BH binaries due to a starburst in the 
early universe. However, as discussed in PS-GRB, the obvious subsets of simulations for which wc can reliably make 
predictions, such as the set of simulations with n > 3000, could have been biased, over-representing those simulations 
which predict more mergers per simulated binary. To mitigate the influence of this bias, in addition to requiring that 
each simulation have many BH-BH binaries (n > 3000), we also (or, equivalently, instead) require that nN > 1 x 10^. 
The second cut's prefactor (1 x 10^) is chosen so that every simulation that satisfies it has n > 3000. The subset 
of simulations which satisfies both conditions has a relatively unbiased distribution of n/N , the fraction of simulated 
binaries that end up as bound BH-BH binaries (PS-GRB). These two cuts do not introduce biases in any of the 
parameters: the set of simulations which satisfy these cuts has similar statistical properties as the entire archive (we 
omit the many plots needed as proof). 

Handling undersampled cases: Finally, with each simulation guaranteed to have 3000 binaries, most simulations have 
a minimum of ^ 1 merging BH-BH binaries; see Figure [TT] for details. Some simulations still do contain very few 
binaries that are merging or even nearly so, largely because certain combinations of physical circumstances rarely lead 
to the formation of many tight black hole binaries. Given that we construct model universes from independent ellptical 
and spiral simulations, these simulations will lead to an indeterminate, small detection rate in only 1/NeNs — 1/300^ 
universe models, multiplied by a factor of order unity. Even considering only elliptical or spiral galaxy star formation, 
they influence our predictions 0{1/Ne,s) — 1% of the time, multiplied by a factor of order unity. We therefore adopt 
a few percent as the lowest resolvable probability of any event our simulations can resolve. 

Potential biases in multiply-selected simulations: Just as the set of simulations was reduced to find simulations with 
many BH-BH binaries in § [H here we require our simulations additionally have enough NS-NS binaries to enable 
accurate estimatesEH Most population synthesis simulations wc performed that produce an adequate number of BH- 
BH binaries also produce more than enough NS-NS and NS-BH binaries 10^) to allow us to estimate the three 
properties {X, dP/dt,p{Aic) needed to reconstruct merger and detection rates for these NS binaries. However, certain 
combinations of population synthesis parameters strongly favor merging BH-BH production over NS-NS production. 
Most of our simulations of these conditions did not produce enough NS-NS binaries to allow accurate estimates. In 
other words, the set of simulations that have both enough NS-NS and close BH-BH binaries is slightly biased against 
the highest BH-BH and lowest NS-NS merger rates; sec Table [T] 

The goal of this paper is to determine the relative likelihood of CBC merger and detection rates given existing 
observations. In particular, in ? l5.2l we examine the distribution of BH-BH merger rates for (spiral-galaxy) simulations 
for which the NS-NS merger rates are consistent with milky Way observations. This subset of simulations is inevitably 
contained within the biased set of simulations for which NS-NS and BH-BH merger rates can be accurately estimatated. 
That being said, the observational constraints reviewed in § 15.21 point to a high NS-NS merger rate and will therefore 
inevitably rule out the few models for which NS-NS merger rates could not be accurately determined. 



We adopt the same conditions on n and TV as were used in 



PS-GRB. 



